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1  Executive  Summary 

The  Pulse  Power  Hybrid  Energy  Storage  Module  Development  Program  compared  two 
systems  to  power  large  pulsed  loads.  The  notional  load  used  in  the  comparison  was  a  high 
energy,  capacitor-based,  Pulse  Forming  Network  (PFN)  operated  with  a  repetition  rate  of  12 
charge/discharge  cycles  per  minute.  Two  hybrid  solutions  were  evaluated.  One  was  based  on 
lithium  ion  batteries  and  used  a  flywheel  to  improve  the  power  delivery.  The  other  was  a 
rotating-machine  based  Hybrid  Energy  Storage  Module  (HESM)  designed  to  mitigate  the  impact 
of  transient  (pulsed)  loads  on  the  ship’s  distribution  power  system.  An  important  finding  was  that 
the  transient  load  could  be  effectively  buffered  with  only  the  flywheel  energy  storage  element. 


Two  PFN  charging  profiles  were  evaluated;  constant  current  and  constant  current/constant 
power.  The  constant  current  profile  was  determined  to  be  the  more  efficient  method  for  charging 
capacitors  bur  requires  a  peak  charging  power  of  twice  the  average  power.  Hybrid  constant 
current/constant  power  charging  is  the  recommended  profile  because  it  provides  high  efficiency 
but  reduces  the  peak  power  demands  on  the  charging  power  supplies  and  energy  storage  sub¬ 
system. 
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Two  topologies  for  the  HESM  placement  within  the  ship’s  power  system  were  considered; 
centralized  placement  of  a  single  HESM  co-located  with  the  PFN  and  distributed  architectures 
with  multiple  HESMs  placed  in  different  zones  of  the  power  distribution  system.  Each  topology 
was  evaluated  using  both  of  the  PFN  charging  profiles.  The  PFN  charging  profiles  are  simulated 
first  using  controlled  ideal  sources  for  the  HESM  and  the  PFN.  The  ideal  results  showed  that  the 
approach  was  sufficiently  appealing  that  more  detailed  assessments  were  warranted.  The 
“actual”  assessments  used  more  detailed  models  for  the  power  electronic  converters,  generators 
and  motors  used  within  the  HESM  and  the  ship’s  power  system. 
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Organization  of  the  HESM  Study 


Control  of  the  HESM  was  accomplished  by  means  of  derived  voltage  reference  signals  that 
served  to  control  the  power  flow  in  and  out  of  the  HESM  controller  and  to  maintain  constant  ship 
current.  The  simulation  results  validated  this  as  an  effective  approach  for  control  of  the  HESM’s 
bi-directional  converter. 

The  flywheel  power  and  energy  requirements  closely  match  the  characteristics  of  an 
advanced  pulsed  power  supply  developed  under  the  Combat  Hybrid  Power  Supply  program 
(CHPS).  A  reassessment  of  this  topology  showed  that  it  should  meet  the  electrical  requirements 
and  the  size  requirements  with  some  modification  of  the  earlier  design  that  was  optimized  for 
somewhat  different  requirements.  Development  and  integration  studies  for  a  customized  Navy 
variant  of  the  concept  --  the  CHPS-N  machine  —  are  being  conducted  under  a  separate  ONR 
grant  (ONR  Grant  N00014-14-1-0262). 

The  work  performed  showed  that  an  architecture  based  on  the  HESM  concept  to  minimize 
the  impact  of  pulsed  power  loads  on  a  ship’s  power  system  is  very  effective.  Furthermore  it 
confirms  that  the  HESM  can  be  used  as  an  additional  support  element  for  the  ship’s  power 
system  when  not  needed  for  pulsed  load  charging  (e.g.  UPS  duty),  thus  making  it  a  versatile 
component  for  a  variety  of  tasks. 
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2  Background 

Future  pulsed  weapon  and  sensors  will  create  challenges  to  Navy  shipboard  power  systems 
[1,  2]  as  transient  loads  will  constitute  a  significant  fraction  of  the  installed  capacity.  The  goal  of 
this  research  program  is  to  explore  the  design  and  integration  of  a  Hybrid  Energy  Storage 
Module  (HESM)  that  can  work  to  mitigate  the  impact  of  transient  (pulsed)  loads  on  the  ship’s 
power  distribution  system. 

The  evaluation  has  been  conducted  in  the  context  of  a  high-energy  capacitor-based  Pulse 
Forming  Network  load  (PFN)  operated  with  a  repetition  rate  of  12  charge/discharge  cycles  per 
minute.  The  PFN  consists  of  multiple  capacitor  modules  served  by  independent  dc  charging 
power  supplies.  Since  this  research  focuses  on  the  HESM,  modeling  costs  were  reduced  by 
representing  the  PFN  as  a  single  aggregated  load. 

The  baseline  PFN  charge/discharge  profile  requires  delivering  80  MJ  of  energy  at  80  %  duty 
cycle  with  a  five  second  repetition  rate.  This  translates  to  a  20  MW  pulsed  load  on  the  ship 
power  system  for  four  seconds  followed  by  a  one  second  break  to  discharge  the  PFN  (1 6  MW 
averaged  over  5  seconds).  (See  green  waveform  on  the  left  side  of  Fig.  4.)  This  duty  cycle  and 
load  level  would  apply  unacceptable  transient  loads  to  the  ship  power  system.  Therefore,  the 
goal  of  the  HESM  is  to  provide  load  leveling  for  the  pulsed  PFN  charging. 

There  are  two  primary  options  for  charging  of  a  capacitive  energy  store:  ramped 
power/constant  current'  and  hybrid  power.^  Constant  current  charging  is  the  more  efficient 
method  far  charging  capacitive  energy  stores  and  is  the  preferred  approach  for  capacitor 
charging  at  low  power  levels.  In  this  scheme,  the  current  into  the  capacitive  store  is  constant  and 
the  power  provided  by  the  charging  power  supply  increases  linearly  with  capacitor  voltage.  This 
results  in  a  ramping  power  profile  with  a  peak  of  approximately  twice  the  average  power.  For 
this  application,  constant  current  charging  requires  a  peak  charging  power  of  ~40  MW  and  a 
charging  system  sized  for  this  load  even  though  it  represents  only  a  small  fraction  of  the  duty 
cycle. 

Simulations  demonstrate  that  the  hybrid  power-charging  scheme  requires  slightly  higher 
currents  from  the  PFN  charging  power  supplies,  but  the  peak  power  demand  on  the  charging 
power  system  is  much  lower.  Furthermore,  in  the  hybrid  approach  proposed  herein,  the  charging 
power  profile  is  initially  ramped  before  it  transitions  to  a  constant  power  mode  that  reduces  the 
peak  demand  on  the  ship  power  system  and  HESM. 


'Ramped  power  upstream  of  the  PFN  power  supply  corresponds  to  constant  capacitive-charging-current  downstream 
of  the  PFN  power  supply. 

^  Power  delivered  to  PFN  power  supply  first  ramps  then  remains  constant. 
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3  Report  Organization 

This  report  is  organized  as  depicted  by  Fig.  1.  The  research  addressed  two  topologies  for  the 
HESM  placement:  centralized  and  distributed.  Subsequently,  each  topology  is  evaluated  using 
two  PEN  charging  variants:  ramped  power  and  hybrid  power.  Each  PEN  charging  profile  is  then 
simulated  using  two  approaches.  When  using  ideal  sources  for  the  HESM  and  PEN,  the 
simulation  type  is  termed  ideal  to  distinguish  them  from  the  case  of  using  a  power  converter  and 
machine  model  to  represent  the  HESM  (termed  actual).  The  assessment  was  initially  performed 
with  idealized  components  as  a  screening  tool.  The  ideal  results  showed  that  the  approach  was 
sufficiently  appealing  that  more  detailed  assessments  were  warranted.  Differences  between  the 
ideal  and  actual  simulations  approaches  are  discussed  as  they  are  presented. 
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Fig.  1.  Organizational  chart  outlining  information  contained  in  this  report.  Level  1:  HESM  placement 
(centralized  vs.  distributed);  Level!:  charging  mode  (ramped  vs.  hybrid);  Level 2^:  simulation  type  (ideal  or 
actual).  Case  studies  presented  herein  are  labeled  (a)  through  (i). 


The  bottom  row  in  Fig.  1  shows  letters  to  indicate  that  this  report  contains  eight  simulation 
result  sets  {a  through  /)  corresponding  to  actual  and  ideal  simulations  of  the  two  HESM 
topologies  and  two  PEN  charging  profiles.  Eor  each  of  these  baseline  case  studies,  a  medium 
voltage  dc  (MVDC)  distribution  system  [3]  was  selected;  alternative  electric  power  distribution 
architectures  would  primarily  impact  the  power  conversion  modules  providing  the  interface  to 
the  HESM  modules. 
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This  remainder  of  this  report  is  organized  as  follows. 

•  Section  4  introduces  the  HESM  topologies 

•  Section  5  introduces  the  PEN  charging  variants 

•  Section  6  presents  simulation  results  of  the  centralized  HESM  topology 

•  Section  7  presents  simulation  results  of  the  distributed  HESM  topology 

•  Section  8  provides  insights  on  the  intended  machine  characteristics 

•  Section  Error!  Reference  source  not  found,  presents  interim  summary  and  conclusions 

•  Section  Error!  Reference  source  not  found,  lists  items  currently  under  investigation 

HESM  Topologies 
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4  HESM  Topologies 


4.1  HESM  Topology  1:  Centralized 

A  one-line  diagram  [4]  of  the  MVDC  distribution  architecture  with  a  centralized  HESM 
topology  is  shown  in  Fig.  2.  The  PEN  and  HESM  are  connected  from  zone  1  at  an  isolated  dc 
bus.  This  topology  implies  a  single  dc  feeder  to  serve  both  the  HESM  and  PEN,  compartmental 
proximity  between  the  HESM  and  PEN,  and  the  maximum  storage  requirement  for  the  single 
HESM  unit.  _ 

HESM  Topology  1:  Centralized 
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Graphic  symbols  for  electrical  and  electronics  diagrams 
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Fig.  2.  JM  VDC  architecture  with  a  centralized  HESM  and  PFN  connected  at  an  isolated  dc  bus. 


A  single  compartment  for 
both  PFN  and  HESM  units 


EHESM  flywheel 


4.2  HESM  Topology  2:  Distributed 

The  one-line  diagram  of  the  MVDC  distribution  architecture  with  a  distributed  HESM 
topology  is  shown  in  Fig.  3.  The  two  HESM  units  considered  for  this  topology  are  distributed 
over  zones  2  and  3.  Additional  HESM  modules  can  be  considered  using  the  same  control 
approach  developed  for  this  initial  case.  The  PFN  remains  connected  to  zone  1  for  consistency 
with  the  centralized  topology  shown  in  Fig.  1.  The  advantages  of  the  distributed  topology  are 
the  inter-zonal  placements  that  allow  the  HESMs  to  supplement  generation  ride-though  when  the 
ring  bus  is  open,  reduced  failure  probability  under  battle  impact,  and  simplified  integration. 
Furthermore,  each  distributed  HESM  unit  is  smaller  and  weighs  less  than  the  single  HESM  unit 
in  the  centralized  case. 
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Fig.  3.  MVDC  architecture  with  distributed  (zonal)  HESMs  and  a  PFN. 


The  assumed  lengths  and  per-unit-length  impedances  for  the  relevant  shipboard  cables  and 
feeders  used  in  the  simulation  models  are  listed  in  Table  1.  As  a  first  approximation,  conductor 
types,  ampacities,  voltage  regulation,  material  characteristics,  flexibility,  geometric  shape, 
economics,  operating  temperature,  insulation  [5],  standing,  temperature,  and  number  shipboard 
cable  conductors  [1,  6,  7]  required  to  supply  the  necessary  current  levels  are  not  considered  as 
they  have  a  minimal  impact  on  the  simulation  results.  However,  conductor  selection  is 
contingent  upon  these  considerations  [8]  which  may  be  important  in  more  detailed  studies. 

Table  1 .  Resistance  of  relevant  cables. 


Cable  < 

or  feeder 

Length 

(feet) 

DC 

Resistance* 

(Cl) 

From 

To 

Ship  ring 

Isolated  dc  bus 

16.4 

133E-6 

Isolated  dc  bus 

HESM 

32.8 

266E-6 

Isolated  dc  bus 

PFN 

32.8 

266E-6 

Ship  ring 

HESM  unit  1 

16.4 

133E-6 

Ship  ring 

HESM  unit  2 

16.4 

133E-6 

Port 

Starboard 

328 

3E-3 

*Calculated  at  8.1 1  mn/1000  ft.  at  25  °C 
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5  PFN  Charging 

In  the  centralized  case  of  Fig.  2,  the  ship  power  {Pship)  is  the  power  delivered  from  the  ship’s 
ring  bus  to  the  dc  bus  connected  to  it  (zone  1).  The  objective  of  the  HESM  is  to  minimize 
transients  on  this  parameter.  The  PFN  power  (Ep/«)  is  the  power  delivered  directly  to  the  PFN 
power  supply.  The  HESM  has  two  power  levels  associated  with  it;  a  positive/charging  power 
and  a  negative/discharging  power.  Both  of  these  power  levels  are  denoted  as  Phesm  and  vary  only 
in  their  direction  (positive  to  mean  charging  or  negative  to  mean  discharging).  The  HESM 
charging  power  (Phesm>0)  is  the  power  delivered  from  the  isolated  dc  bus  to  the  HESM  bi¬ 
directional  converter  as  shown  in  Fig.  2.  This  power  flows  through  the  bi-directional  converter 
and  motor/generator  toward  the  flywheel.  The  HESM  discharging  power  (Phesm<0)  is  the  power 
flowing  away  from  the  HESM  bi-directional  converter  toward  the  isolated  dc  bus.  The  HESM 
discharge  power  supplements  the  ship  power  when  the  PFN  charges — that  is,  when  the  PFN 
charges  Ppf,i=Pship-Phesm\  and  when  the  PFN  idles  and  the  HESM  charges,  Pp/„=0  and  Pship^Phesm- 
Note  that  Phesm  is  negative  when  the  HESM  is  supplying  power  to  the  PFN. 

In  the  distributed  case  of  Fig.  3,  Pship  is  the  power  contribution  of  the  two  ship  generators 
toward  charging  the  HESM  and/or  PFN.  The  PFN  power  is  the  power  delivered  from  the  ship 
ring  bus  to  the  PFN  power  supply.  The  total  HESM  power  (Phesm)  is  the  power  delivered  to  or 
from  the  ring  bus  to  the  two  HESM  units  together.^  That  is,  the  total  HESM  power  is 
P hesm=P hesmX'^P hesmi-  Similar  to  as  shown  in  Fig.  2,  Phesm  is  shown  in  the  charging  direction  (into 
the  HESM). 

5.1  PFN  Charging  Variant  1:  Ramped  Power/Constant  Current 

The  first  PFN  charging  variant  considered  is  the  ramped  power/constant  current.  The  PFN 
charging  system  can  be  modeled  as  an  array  of  independent  high-voltage  charging  power 
supplies  each  charging  individual  capacitor  modules.  This  charging  scheme  is  the  most  efficient 
but  requires  the  highest  peak  power  from  the  PFN  charging  power  supplies  and  HESM.  The 
baseline  ship,  HESM,  and  PFN  power  profiles  corresponding  to  the  ramped  power  charging 
profile  are  shown  on  the  left  of  Fig.  4.  The  goal  of  the  HESM  is  to  level  the  PFN  charging  load 
seen  by  the  ship  power  system — that  is,  to  level  Pship  to  a  near  constant  value.  During  the  initial 
stages  of  the  PFN  charging  profile  where  power  into  the  PFN  is  low,  the  HESM  accepts  power 
from  the  ship  (Phesm  >  0)  indicating  charging  of  the  energy  store.  As  the  PFN  power  increases 
and  reaches  Ppf„=Pship,  the  HESM  power  reverses  to  provide  power  (Phesm^O)  to  the  isolated  dc 
bus  to  supplement  Pship  during  the  latter  stages  of  the  PFN  charging  profile,  which  is  when  Ppfn  > 

P ship- 

5.2  PFN  Charging  Variant  2:  Hybrid  Power 

The  second  PFN  charging  variant  is  a  hybrid  power  scheme  that  results  in  the  power  profiles 
shown  on  the  right  of  Fig.  4.  The  PFN  charges  with  a  combination  of  constant  capacitive- 
charging  current  (rising  power  edge  between  t  =  0  and  t=\  s)  and,  subsequently,  constant  power 
(between  t  =  1  and  t  =  4  s). 


^  It  is  assumed  that  the  HESM  units  charge  and  discharge  in  unison. 
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Under  this  second  charging  variant,  the  discharge  HESM  power  is  different.  In  charging 
variant  1,  the  HESM  discharge  power  was  Phesm~^^  MW  while  in  variant  2  it  is  Phesm-^  MW. 

The  charging  HESM  power,  however,  in  both  variants  is  the  same.  The  ship  power  in  both  cases 
is  also  the  same  for  variants  1  and  2,  which  satisfies  the  aforementioned  goal  of  ship  load¬ 
leveling  with  an  HESM.  The  primary  advantages  of  charging  variant  2  are  reduced  peak  HESM 
power,  reduced  energy  transfer  through  the  HESM,  and  reduced  peak  power  from  the  charging 
power  supplies. 
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Power  (MW) 


PFN  Charging  Variant  1: 
Ramped  Power  (Constant  Current) 


PFN  Charging  Variant  2: 
Hybrid  Power 

Pship^llMW 


Fig.  4.  Intended  power  profiles  for  the  ship,  PFN,  and  HESM  for  a  centralized  topology.  Left:  ramped  PFN  charging.  Right:  hybrid  charging.  The  HESM  power  (blue) 

levels  the  ship  power  demand  to  17  MW  (red)  during  the  PFN  charging  cycles  (green). 
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6  Centralized  HESM 

This  section  analyses  the  performance  of  the  HESM  unit  in  the  centralized  topology  shown  in 
Fig.  2.  A  voltage  reference  signal  for  the  HESM  bi-directional  converter  is  derived  first.  This 
voltage  reference  signal  allows  the  system  to  charge  and  discharge  the  flywheel  energy  store 
while  maintaining  the  ship  power  constant  (leveled).  The  voltage  reference  signal  derived  next 
requires  specifying  the  desired  (leveled)  ship  current,  and  will  be  derived  using  an  amp-hour 
balance.  (Alternative  control  schemes  are  possible;  this  approach  was  selected  to  allow  the  study 
to  focus  on  characterizing  the  rotating  machines.) 


After  the  derivation  of  the  voltage  control  signal,  the  simulation  results  validate  the  HESM 
bi-directional  converter  voltage  controller  against  the  intended  profiles  shown  on  the  left  of  Fig. 
4.  The  simulation  results  include  the  energy  stored  (delivered)  to  (from)  the  flywheel,  which 
brings  about  important  considerations  to  size  the  HESM  motor/generator. 

6. 1  Centralized  Voltage  Control 


This  section  derives  a  voltage  reference  signal  for  the  HESM  bi-directional  converter. 
Referring  to  Fig.  2,  the  sum  of  currents  at  the  isolated  dc  bus  is  given  by  (6. 1),  where  fhip  is  the 
ship  current  flowing  from  the  ring  bus  to  the  isolated  dc  bus,  fesm  is  the  HESM  charging  current, 
and  Ipfn  is  the  PEN  charging  current.  As  it  is  noted  from  Fig.  2  and  Fig.  3,  these  currents  follow 
the  directions  of  their  respective  powers. 


^ ship  ^hesm~^^pfn  (^-l) 

Expressing  fesm  in  terms  of  the  isolated  dc  bus  voltage  Vdc  and  HESM  converter  output  voltage 
Vhesm  results  in  (6.2),  where  Rcabie  is  cable  (feeder)  resistance  between  the  ring  and  isolated  dc 
bus. 


...  {V.c-V„esn,) 


Substituting  (6.2)  into  (6.1)  results  in  (6.3),  where  solving  for  the  HESM  converter  output 
voltage  results  in  (6.4). 


fhip 


{V.c -Vt.es.) 


R 


+  / 


pfn 


cable 


V  ^  V  - 

hesm 

must  be 
measured 


R  (I 

cable  ship 


must  be 
measured  j 


(6.2) 


(6.3) 

(6.4) 


Equation  (6.4)  is  the  voltage  reference  signal  of  the  HESM  bi-directional  converter.  This 
reference  signal  determines  whether  the  HESM  should  charge  (when  Vdc>  Vhesm  — ^  Phesm>0)  or 
discharge  (when  Vdc<Vhesm  Phesm<^)-  This  equation  also  indicates  that  Vdc  and  Ipf„  must  be 
measurable  quantities  and  that  fkip  is  an  exogenous  constant,  which  is  derived  using  an  amp-hour 
balanced  in  section  6.2.1 . 

Specifying  fkip  requires  exercising  care  as  inaccurate  values  result  in  the  HESM  flywheel 
drifting  in  speed.  For  example,  if  the  HESM  charging  amp-hours  are  positively  unbalanced,  the 
HESM  flywheel  will  continue  to  (slowly)  accelerate  indefinitely.  If  the  HESM  charging  amp- 
hours  are  negatively  unbalanced,  the  HESM  flywheel  will  (slowly)  decelerate  until  it  stalls. 
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6.2  Ramped  Power 

This  section  first  derives  a  value  for  Iship  in  (6.4)  using  an  amp-hour  balance.  Following,  the 
simulation  results  for  the  HESM  centralized  topology  with  ramped  power  /  constant  current  PFN 
charging  are  presented  (result  set  (a)  in  Fig.  1). 

6.2.1  Amp-Hour  Balance 

Referring  to  the  HESM  charging  power  on  the  left  of  Fig.  4  (blue  waveform),  the 
corresponding  charging  current  Ihesm  over  one  charging  cycle  is  depicted  with  Fig.  5  .  When 
hesm^Q,  the  HESM  charges;  when  lhesm<0,  the  HESM  discharges.  This  profile  shows  three  time 
segments  denoted  by  T],  T2,  and  T3.  During  Ti,  Ihesm  (initially  charging  the  HESM  at  hesm^hhip) 
decays  linearly  from  to  .  The  ending  value  of  hesm  is  negative, 

which  means  that  Ihesm  reverses  direction  to  discharge  the  HESM  when  Ipfn>Iship-  During  T2,  hesm 
reduces,  reverses,  and  ramps  back  to  the  maximum  HESM  charging  current  of  hesm—Iship-  During 
Tj,  the  HESM  current  is  maintained  at  hesm-hhip  to  re-charge  the  HESM.  This  three-segment 
cycle  repeats  every  5  seconds. 


Fig.  5.  Amp-hour  balance  for  ramped  power  (charging  variant  1). 


Referring  to  Fig.  5,  the  desired  value  of  Imax^hhip  is  approximated  by  balancing  the  shaded 
amp-hour  areas.  The  net  sum  of  amp-hours  is  given  in  (6.5),  where  solving  for  hhip  results  in 
(6.6). 


1 


1 


ship 


-I 


-t  +  I 


dt- 


(6.5) 


0.425 


(6.6) 


Equation  (6.6)  suggests  that  to  maintain  constant  ship  power  demand  and  prevent  flywheel  speed 
drifts,  Iship  should  be  ~40%  of  the  peak  PFN  charging  current.  Referring  to  the  time  periods 
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shown  in  Fig.  2,  substituting  Ti=4  s,  T2=0.25  s,  73=0.75  s,  and  =  40  MW/6  kV  «  6,666  A 
in  (6.6)  results  in  (6.7). 

7,„«  2,833  A  (6.7) 

The  constant  (leveled)  value  of  Iship  in  (6.7)  must  be  maintained  by  the  ship  to  prevent 
flywheel  speed  drifts.  This  value  conforms  to  the  expected  ship  power  demand  of  Pship  ~  17  MW 
on  the  left  of  Fig.  4  and  as  given  by  (6.8). 

=  2.83  kA  X  6  kV  « 1 7  MW  (6.8) 

Having  derived  Ighip  for  use  with  the  voltage  reference  signal  of  (6.4),  the  following 
subsection  presents  preliminary  simulations  results  to  validate  (6.4)  and  (6.7). 

6.2.2  Simulation  Results 

This  subsection  presents  simulation  results  of  the  scenario  in  Fig.  2  using,  first,  power- 
controlled  ideal  sources  for  the  PFN  and  HESM  units  ((a)  and  (c)  in  Fig.  1)  and,  then,  using  a  bi¬ 
directional  converter,  motor/generator,  and  flywheel  for  the  HESM  ((b)  and  (d)  in  Fig.  1).  These 
two  simulation  run  sets  are  termed  ideal  and  actual  in  Fig.  1,  respectively. 

6.2.2. 1  Simulation  with  Ideal  Sources  (a) 

Referring  to  Fig.  2,  the  ring  bus  was  modeled  with  two  rectified  synchronous  generators 
serving  10  MW  zonal  loads.  The  PFN  was  modeled  as  a  power-controlled  current  source 
calculated  during  run  time  as  I =  PpfnjVpp^ ,  where  Ppfa  is  known  from  the  power  profile  on  the 

left  of  Fig.  4  and  Vpf„  is  the  voltage  measured  at  the  PFN  input  terminals.  The  HESM  was 
modeled  as  a  controlled  voltage  source  driven  with  (6.4)  and  (6.7).  The  ideal  simulation  results 
for  the  centralized  HESM  topology,  ramped  power  PFN  charging  profile  are  shown  on  the  left  of 
Fig.  7  and  summarized  next: 

Voltages:  The  top  curve  set  shows  dc  bus,  PFN,  and  HESM  terminal  voltages.  These 
voltages  are  nearly  the  same  as  the  voltage  drop  across  cables  of  short  length  is  negligible. 
Although  not  readily  clearly  visible,  the  HESM  bi-directional  converter  controller  ramps 
Vhesm  about  Vdc  Using  (6.4)  to  produce  bidirectional  power  flows  in  and  out  of  the  HESM. 

The  PFN  voltage  Vpf„  ramps  down  during  the  PFN  charging  periods  due  to  the  voltage  drop 
of  its  feeder,  and  it  returns  to  Vdc  when  the  PFN  idles  between  discharge  cycles. 

Current:  The  second  curve  set  on  the  left  of  Fig.  7  shows  the  PFN,  HESM,  and  ship  current. 
The  leveled  ship  current  is  consistent  with  the  computed  value  of  the  amp-hour  balance  in 
(6.7). 

Power:  The  third  curve  set  on  the  left  of  Fig.  7  shows  Pship  ~17  MW,  which  is  consistent  with 
the  calculation  in  (6.8)  and  with  the  leveled  ship  power  on  the  left  of  Fig.  4  (red  trace).  The 
PFN  charging  power  Ppf„  is  also  in  agreement  with  charging  variant  1  anticipated  by  Fig.  4. 
The  HESM  power  Phesm  is  bi-directional:  Phesm>^  indicates  the  HESM  is  charging  and 
PhestnFt  indicates  the  HESM  is  discharging. 

Energy:  The  fourth  curve  set  on  the  left  of  Fig.  7  shows  the  energy  delivered  from  the  ship 
ring  bus  and  the  energy  consumed  by  the  PFN.  The  energy  stored  in  the  HESM  is  displayed 
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as  flywheel  speed  (RPM)  against  the  right-side  vertical  axis.  The  flywheel  speed  was 
calculated  from  the  relation  of  energy  stored  in  a  rotating  mass  as 

W  =  ]^Jco^  (6.9) 


where 

fV  =  energy  stored  in  flywheel  mass  [Joules] 

J  -  flywheel  moment  of  inertia  [kg-m^] 
o)  =  flywheel  angular  speed  [rad/s]. 

For  the  initial  simulations,  the  flywheel’s  moment  of  inertia  was  set  to  J=l,128  kg-m^,  which 
is  the  same  as  the  rotor  inertia  of  an  Electromagnetic  Aircraft  Launch  System  (EMALS)  pulse 
power  generator  [9].  Expressing  (o  as  revolutions  per  minute  (N,  RPM)  results  in  (6.10). 

W  =  (RPM)  (6.10) 

As  noticed  from  the  HESM  speed  trace  (Nhesm)  in  Fig.  7,  there  is  a  drift  in  flywheel  speed. 
This  drift  is  due  to  a  minor  unbalance  neglected  in  the  amp-hour  analysis.  The  balancing 
equations  do  not  account  for  some  real  life  issues  like  dc  voltage  ripple,  the  exact  voltage  level  at 
the  isolated  dc  bus,  and  the  exact  charging/discharging  values  of  Ihesm-  The  flywheel  speed  drift 
has  been  in  fact  compensated  by  the  addition  of  a  speed  control  loop  feedback  on  the  power 
converter.  This  speed  loop  feedback  was  added  in  the  final  version  of  the  models  delivered  to  the 
Navy  and  has  reduced  the  flywheel  speed  drift  to  an  inconsequential  level.  This  solution,  of 
course,  was  accomplished  by  intervening  via  the  software  on  the  model  design  not  on  the 
foundational  analytical  equations  described  above  which,  in  and  of  themselves,  do  not  address 
this  issue. 


6.2.2.2  Simulation  with  Actual  Sources  (b) 

Referring  to  Fig.  2,  the  two  generators  providing  the  power  to  the  main  ring  bus  were 
modeled  as  wound  field  synchronous  machines  with  their  output  rectified  by  a  12-pulse  rectifier. 
The  field  voltage  of  the  synchronous  generator  was  regulated  to  maintain  a  constant  output 
voltage  of  4,444  VAC  (rms,  line-to-line)  rectified  to  6  kV  dc. 

The  HESM  motor/generator  was  modeled  as  a  programmable  three-phase  voltage  source 
with  a  nominal  output  of  3,500  V  ac  (rms  line-to-line)  and  with  an  inertia  of  1,128  kg-m^ 
(corresponding  to  the  EMALS  machine)  representing  the  flywheel  energy  storage.  The  flywheel 
storage  can  be  either  the  rotor  of  the  motor/generator  or  a  separate  rotating  unit  coupled 
mechanically  to  it.  The  motor/generator  is  connected  to  the  isolated  dc  bus  via  a  bi-directional 
power  converter.  The  converter  accepts  variable  frequency  ac  power  from  the  flywheel- 
generator  storage  and  provides  rectified  output  dc  power,  or  accepts  power  from  the  isolated  dc 
bus  and  converts  it  to  three-phase  ac  power  to  drive  the  synchronous  machine  as  a  motor  and 
store  energy  in  the  flywheel.  The  bi-directional  converter  is  controlled  by  considering  the  d-q 
axes  representation  of  its  current  and  actively  controlling  the  direct  axis  component  [10].  The 
PFN  load  was  modeled  as  a  power-controlled  current  source  just  as  was  done  in  the  ideal  case. 

The  results  of  the  simulations  with  actual  sources  are  summarized  in  Fig.  9  (left).  As  can  be 
seen,  they  reproduce  very  closely  the  results  obtained  with  ideal  sources  (Fig.  7)  except  for  the 
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voltages.  In  the  ideal  simulations  ((a)  and  (c)),  the  six-pulse  generator  rectifiers  and  absence  of 
voltage  regulation  produced  significant  dc  voltage  ripple.  However,  this  ripple  was  ignored 
because  it  did  not  affect  the  simulation  results  and  because  it  is  corrected  in  the  actual  simulation 
result  set  ((b)  and  (d)).  The  following  observations  can  be  made  about  Fig.  9: 

Voltage:  The  PFN  voltage  Vpf„  in  Fig.  9  (left)  practically  coincides  with  the  ship’s  bus 
voltage  Vdc  (green  trace  is  beneath  red  trace)  while  the  HESM  bi-directional  converter 
controller  ramps  Vhesm  about  Vdc  according  to  (6.4)  to  produce  bidirectional  power  flows  in 
and  out  of  the  HESM. 

Current:  The  ship’s  current  is  kept  nearly  constant  during  the  PFN  charge-discharge  cycle 
although  it  varies  non-linearly  between  a  minimum  of  2,550  A  and  a  maximum  of  2,950  A 
around  an  average  of  about  2,850  A. 

Power:  The  PFN  and  HESM  traces  here  mirror  those  obtained  in  the  ideal  case.  The  ship’s 
power,  however,  is  not  exactly  constant  but  varies  non-linearly  between  a  minimum  of  15.5 
MW  and  a  maximum  of  1 8  MW  around  an  average  of  about  17  MW  corresponding  to  the 
current  variation. 

Energy:  The  traces  here  are  very  similar  to  those  obtained  in  the  ideal  case.  The  only 
observation  relates  to  the  HESM  stored  energy  (flywheel  RPM)  which  here  drifts  downward 
whereas  in  the  case  of  ideal  sources  it  tended  to  increase.  As  was  observed  previously,  this 
drift  is  due  to  minor  unbalances  neglected  in  the  amp-hour  analysis  and  can  be  eliminated  by 
proper  tracking  of  N  and  correction  of  Iship  in  (6.4). 

6.3  Hybrid  Power 

This  section  presents  simulation  results  for  the  HESM  centralized  topology  using  the  hybrid 
power  PFN  charging  shown  on  the  right  of  Fig,  4.  Similar  to  the  constant  power  case,  before 
presenting  the  simulation  results,  a  value  for  Chip  in  (6.4)  must  is  derived  according  to  the 
following  an  amp-hour  balance. 

6.3.1  Amp-Hour  Balance 

Because  the  hybrid  power  variant  considers  the  same  centralized  HESM  topology  of  Fig.  2, 
the  HESM  voltage  control  signal  derived  in  (6.4)  remains  the  same  as  in  the  centralized  case. 

The  value  for  the  ship  current  Chip  in  (6.4),  however,  must  be  re-calculated  as  the  waveform  of 
hesm  (Fig,  6)  in  the  hybrid  charging  profile  is  different  than  it  was  for  the  ramped  power  charging 
profile  (Fig,  5). 

The  HESM  charging-current  profile  Cesm  over  one  charging  cycle  of  hybrid  PFN  charging  is 
depicted  in  Fig,  6.  This  profile  shows  four  time  zones  denoted  Tj,  T2,  and  T3.  During  T\,  Cesm 
reduces  from  its  peak  discharge  value  of  /,„j„  =  Chip  ~  >  reverses  direction,  and  ramps  linearly 

to  its  maximum  charging  value  of  =  Chip  •  During  Tj,  Cesm  remains  at  Imax  to  charge  the 

HESM.  During  C,  Cesm  decays,  reverses,  and  ramps  towards  Cm  to  discharge  the  HESM. 

During  Tij,  Cesm  remains  at  Imin  to  supplement  Chip  during  PFN  charging. 
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Fig.  6.  Amp-hour  balance  for  hybrid  power  (charging  variant  2). 


Referring  to  Fig.  6,  the  value  of  Iship  is  estimated  by  balancing  the  shaded  amp-hour  areas  as 
given  by  (6.1 1)-(6.12). 


A  I  -I  ■ 

I  max _ min  *  J 

J  r 


dt  +  f  )  dt  +  j  + 1..,., 


ship 


J  _  jpeak 
^  ship  ^  pfn 


/;;r(-7;-27;-7;) 


0.725  x/T* 


(6.12) 


2(7;  +7;  +7;  +7;) 

Substituting  ri=0.25  s,  72=0.75  s,  73=1  s,  74=3  s,  and  =  22.86  MW/6  kV  «  3,810  A  in 
(6.12)  results  in  (6.13). 

2,762  A  (6.13) 


*  ship 


Iship  in  (6.13)  is  less  than  it  was  for  the  ramped  power  variant  in  (6.7).  Ihe  reduction  in  IsMp 
reduces  the  power  demand  from  17  MW  (eq.  (6.8))  to  16.5  MW  as  given  in  (6.14).  In  both 
cases,  however,  the  ship  power  demand  is  near  the  expected  17  MW  anticipated  by  Fig.  4 

7^,„^=2.762kAx6kV«16.5MW  (6.14) 


6.3.2  Simulation  Results 

Referring  to  Fig.  1,  section  6.2.2  presented  the  result  sets  for  cases  (a)  and  (c).  This  section 
presents  the  result  sets  for  cases  (b)  and  (d). 

6.3.2.1  Simulation  with  Ideal  Sources  (c) 

The  simulation  results  of  the  hybrid  power  PFN  charging  are  shown  on  the  right  of  Fig.  7 
and  are  explained  next. 
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Voltages'.  Comparing  the  left  and  right  side  waveforms  in  Fig.  7,  there  are  no  noticeable 
changes  in  the  voltages  profiles. 

Current'.  The  ship  current  is  similar  in  both  cases;  Iship~  2.7  kA.  The  peak  PFN  charging 
current,  however,  reduced  from  ~6.6  kA  ~3.8  kA.  The  HESM  charging  current  is  the  same 
in  both  charging  variants.  The  peak  HESM  discharge  current,  however,  also  significantly 
reduced  from  —3.8  kA  to  —1.1  kA.  Reductions  in  current  imply  less  interrupting 
requirements  and  cable  sizes,  but  also  sustained  currents  for  longer  times. 

Power:  The  third  curve  set  shows  the  same  approximate  ship  power  demand  as  in  variant  1 : 
Pship  ~  17  MW  (really,  16.5  MW  in  the  hybrid  power  variant).  The  peak  PFN  charging 
power  reduced  from  Ppf„  —  40  MW  to  Ppfr~23  MW,  and  the  HESM  peak  power  reduced  from 
Phesm  ~  23  MW  to  Phesm  ~  17  MW.  A  Contrast  of  the  peak  powers  in  each  PFN  charging 
variant  is  summarized  below: 

•  Ramped  power:  P^hip  ~  17  MW,  Ppf„  -  40  MW,  and  Phesm-^S  MW  (when  discharging) 

•  Constant  power:  P^hip  -  1 6.5  MW,  Ppf„  ~  23  MW,  and  Phesm-^  7  MW  (when  charging) 

Energy:  The  fourth  curve  set  shows  the  energy  delivered  from  the  ship  ring  bus  as  Wship  and 
the  energy  consumed  by  the  PFN.  The  energy  stored  in  the  HESM  is  shown  in  Joules  on  the 
left-side  vertical  axis  and  in  tenus  of  flywheel  speed  (RPM)  on  the  right-side  vertical  axis.  A 
noticeable  difference  in  the  charging  variants  is  seen  in  the  energy  transfer  through  the 
HESM.  Referring  to  Fig.  8,  the  energy  transfer  in  the  hybrid  charging  profile  was  reduced  by 
—29  %  (100x(27. 88-19. 68)727.88=29. 41).  The  hybrid  PFN  charging  (variant  2)  reduces  the 
HESM  peak  power,  energy  transfer  requirement,  and  the  required  HESM  energy  storage. 

6.3.2.2  Simulation  with  Actual  Sources  (d) 

The  results  of  the  simulations  with  actual  sources  are  summarized  in  Fig.  9  (right).  As  can  be 
seen,  they  reproduce  very  closely  the  results  obtained  with  ideal  sources  in  (Fig.  7)  and  the 
following  observations  can  be  made  in  a  manner  similar  to  what  was  done  in  the  constant  current 
charging  case. 

Voltage:  The  PFN  voltage  Vpp  again  practically  coincides  with  the  ship’s  bus  voltage  Vdc 
(green  trace  is  under  the  red  trace)  while  the  HESM  bi-directional  converter  controller  varies 
Vhesm  about  Vdc  according  to  (6.4)  to  produce  bidirectional  power  flows  in  and  out  of  the 
HESM.  The  cycle  is  different  from  the  ramped  power/constant  current,  but  the  isolated  dc 
bus  voltage  remains  constant  around  6  kV  dc. 

Current:  The  ship’s  current  is  kept  nearly  constant  (variations  less  than  ±1%  of  mean)  during 
the  PFN  charge-discharge  cycle  although  it  varies  non-linearly  between  a  minimum  of  2,570 
A  and  a  maximum  of  2,870  A  around  an  average  of  about  2,750  A. 

Power:  The  PFN  and  HESM  traces  here  mirror  those  obtained  in  the  ideal  source  case.  The 
ship’s  power,  however,  is  not  exactly  constant  but  varies  non-linearly  between  a  minimum  of 
15.3  MW  and  a  maximum  of  17.3  MW  around  an  average  of  about  16.4  MW  corresponding 
to  the  current  variation. 
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Energy.  The  traces  here  are  very  similar  to  those  obtained  with  ideal  sources.  The  only 
observation  relates  to  the  HESM  stored  energy  (flywheel  RPM)  which  here  drifts  downward 
whereas  in  the  case  of  ideal  sources  it  tended  to  increase.  This  drift  is  again  due  to  minor 
unbalances  neglected  in  the  amp-hour  analysis  and  can  be  eliminated  by  proper  tracking  and 
correction  (when  necessary). 
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Fig.  7.  Ideal  simulation  results  for  ramped  and  hybrid  PFN  charging  variants  (centralized  HESM).  (Actual  simulation  results  are  shown  in  Fig.  9.) 
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Fig.  8.  Ideal  flywheel  speed  and  energy  for  ramped  and  hybrid  PFN  charging  variants.  (Close-up  of  the  lower  row  of  traces  in  Fig.  7.) 
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(b)  (d) 

•  Centralized  HESM  •  Centralized  HESM 

•  Ramped  PFN  Charging  (variant  1)  •  Hybrid  PEN  Charging  (variant  2) 

•  >Acfaa/ Simulation  Results  j  •  /Ac^ua/ Simulation  Results 


Fig.  9.  Actual  simulation  results  for  ramped  and  hybrid  PFN  charging  variants  (centralized  HESM).  (Ideal  simulation  results  were  shown  in  Fig.  7.) 
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7  Distributed  HESMs 

This  section  addresses  the  right  side  of  the  organization  diagram  presented  in  Fig.  1. 
Following  the  proeedure  illustrated  for  the  centralized  HESM  unit  in  Seetion  6,  voltage  reference 
signals  are  derived  first  for  each  HESM  bi-directional  converter.  After  these  derivations, 
simulation  results  are  presented  to  contrast  the  performance  of  distributed  HESM  units  over  the 
centralized  one.  The  simulation  results  include  the  speed  and  energy  characteristics  of  the 
HESM  flywheels  as  well. 

7. 1  Distributed  Voltage  Control 

Referring  to  the  one-line  diagram  with  distributed  HESM  units  in  Fig.  3,  and  for  sake  of 
analysis,  it  is  assumed  that  the  voltage  around  the  ring  bus  is  the  same.  This  assumption  allows 
representing  the  ring  bus  as  a  single  bus  bar  and,  therefore,  to  simplify  the  one-line  diagram  of 
Fig.  3  to  the  equivalent  in  Fig.  10. 

The  simplified  diagram  in  Fig.  10  preserves  important  eharacteristics  of  Fig.  3.  That  is,  two 
generators  inject  power  to  the  ring  bus,  one  PFN  unit  connects  from  the  ring  bus  on  its  own,  and 
the  two  HESMs  are  regarded  as  distributed  instead  of  centralized.  The  intermediate  bus  between 
the  ship  ring  bus  and  the  HESM  units  is  there  to  show  how  the  total  HESM  power  and  current 
flow  into  each  HESM  unit.  The  simplified  depietion  in  Fig.  10  is  a  eommon  approach  [1 1]  to 
derive  control  expressions  for  each  HESM  bi-directional  converter;  however,  this  simplified 
diagram  is  not  valid  for  simulation.  The  simulations  presented  later  in  this  section  are  based  on 
the  one-line  diagram  shown  earlier  in  Fig.  3. 


current 

Fig.  10.  Simplified  representation  of  ring  bus,  PFN,  and  HESM  units. 

Referring  to  Fig.  10,  it  is  important  to  note  that  Pshtp  and  4/,/p  represent  only  the  contribution 
of  the  ship  generators  toward  the  PFN  and  HESM  units — that  is,  PsMp  and  IsMp  are  not  the  total 
power  and  current  supplied  by  the  ship  generators. 
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The  sum  of  currents  at  the  single-bus  bar  representation  of  the  ring-bus  in  Fig.  10  is  given  by 
(7.1). 

^ ship  ^hesm  ^ pfn  (7.1) 

The  total  HESM  current  (charging  or  discharging)  is: 

T  —J  4-J  —  ^ hem 

hesm  ^hesm\  '^^hesm2  2  2  * 

Substituting  (7.2)  in  (7.1)  results  in: 


'L,I2=!, 


hesrnl  pfhj'^ 
hesm2 


(7.3) 


Referring  to  Fig.  10,  and  following  the  expression  stated  earlier  in  (6.2),  the  current  contribution 
of  each  HESM  can  be  expressed  as  (7.4),  where  Vdc  is  the  voltage  upstream  of  each  HESM  cable, 
Rcabki  is  the  resistance  (Q)  of  the  cable  between  HESM  unit  i  and  its  upstream  interface  bus,  and 
Vhesmi  is  the  desired  (reference)  dc-side  output  voltage  of  HESM  bi-directional  converter  /. 


V 

dc  hesm  2 


) 


,  _ _  ,  _ _  (7.4) 

hesm  1  jp  5  hesm  2  p  ^  ' 

^coble\  ^cablel 

Substituting  (7.4)  in  (7.3)  and  solving  the  voltage  reference  signal  for  each  HEMS  controller 
results  in  (7.6)  and  (7.7)  for  HEMS  units  1  and  2,  respectively. 

{y,c-Vh.s..) 


^  ship 


V  =  V  - 

hesrnl 

must  be 
measured 


n  '  ^  Pfi^ 


-R 


cablel 


V  =  V  -R 

^hesrnl  ^  dc  ^\ahle2 

must  be 
measured 


^  ship 
2 


‘  ship 
2 


must  be 
measured  J 


must  be 
measured  J 


(7.5) 

(7.6) 

(7.7) 


Equations  (7.6)  and  (7.7)  are  the  voltage  reference  signals  of  HESMs  1  and  2,  respectively. 
These  equations  also  indicate  that  each  HESM  will  provide  (share)  one-half  of  lpf„  and  that  each 
HESM  will  charge  at  one-half  of  4/,/p.  This  1/M  current-sharing  technique  [12]  (where  Mis  the 
number  of  HESM  units)  stems  from  the  one -half  HESM  coefficients  introduced  in  (7.2),  which 
can  be  adjusted  as  more  HESM  units  are  considered. 


Similar  to  the  centralized  topology,  the  voltage  controller  of  each  HESM  unit  must  measure 
the  voltage  V^c  at  its  upstream  bus  cormection  and  the  total  PEN  charging  current  lpf„.  The  values 
of  Iship  required  in  (7.6)  and  (7.7)  for  the  ramped  and  hybrid  PEN  charging  power  profiles  were 
derived  previously  with  an  amp-hour  balance  in  (6.7)  and  (6.13),  respectively.  These  values  of 
Iship  are  the  same  in  the  distributed  topology  as  the  energy  delivered  from  the  ship  generators  to 
the  HESMs  and  PEN  does  not  change  with  HESM  topology.  Since  new  values  for  hhip  are  not 
required,  the  simulation  results  are  presented  directly. 
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7.2  Simulation  Results 

Referring  to  Fig.  1,  section  6.2.2  presented  the  result  sets  for  cases  (a)  and  (c).  Section  6.3.2 
presents  the  result  sets  for  cases  (b)  and  (d).  This  section  presents  all  simulation  results  on  the 
right  half  of  the  hierarchical  chart  in  Fig.  1 — that  is,  to  (e),  (f),  (g),  and  (h),  respectively. 

7.2.1  Simulation  with  Ideal  Sources  (e) 

The  ideal  simulation  results  for  the  distributed  HESM  topology  are  shown  in  Fig.  11.  These 
results  are  compared  in  two  ways: 

•  Distributed  vs.  centralized  HESM  (left  of  Fig.  11  vs.  left  Fig.  7,  or  (e)  vs.  (a)) 

•  The  two  PEN  charging  profiles  for  the  distributed  units  (left  vs.  right  of  Fig.  11,  or  (e) 
vs.  (g)) 

7.2.1. 1  Comparison  to  Centralized  Topology 

The  following  explanations  compare  the  left  side  of  Fig.  7  (centralized  HESM)  to  the  left  side 
of  Fig.  11  (distributed  HESMs). 

Voltages:  The  voltage  profile  of  HESM  1  in  Fig.  11  is  closer  to  the  ring  bus  voltage  than  was 
the  centralized  HESM  unit.  This  is  because  there  is  no  longer  a  single  feeder  serving  the 
entire  ship  current  required  by  the  PEN  and  HESM  units.  Instead,  the  feeder  to  each  HESM 
unit  transmits  only  ‘A  of  fhip-  This  lower  current  reduces  the  voltage  drop  in  each  HESM 
feeder. 

Current:  The  magnitude  of  the  ship  current  is  unchanged;  what  changes  is  its  interpretation. 
In  Fig.  7,  Iship  represented  the  ship  current  flowing  in  the  feeder  from  zone  1  to  the  isolated  dc 
bus.  (See  single  feeder  to  isolated  dc  bus  in  Fig.  2.)  In  the  distributed  case  of  in  Fig.  11, 
Iship  represents  the  combined  current  contribution  of  both  generators,  which  serves  each 
HESM  unit  and  the  PEN.  That  is,  fhip  is  not  the  total  current  served  by  each  generator,  fhip 
is  the  current  demanded  by  the  HESM  and  PEN  units  together:  the  current  demanded  by 
other  loads  on  the  ring  bus  will  continue  to  be  provided  by  the  ship’s  generators.  Here  we  are 
focusing  only  on  the  fraction  of  ship’s  current  required  by  the  HESM  and  PEN.  The  PEN 
current,  however,  is  unchanged  as  in  both  centralized  and  distributed  cases  the  PEN  remains 
as  a  single  unit.  Both  the  charging  and  discharging  currents  of  each  HESM  unit  are  now  one 
half  of  the  values  observed  for  the  centralized  HESM  unit.  (In  Fig.  11  all  waveforms  for 
HESM  1  and  2  are  very  similar;  thus,  the  notation  HESM  1 ,2  is  used  to  refer  to  either  HESM 
unit.)  The  reduction  in  HESM  current  demand  favorably  impacts  cable  losses,  voltage  drop, 
ampacity  requirement,  protective  device  sizing,  coordination,  and  placement.  These  impacts, 
however,  are  not  part  of  this  report. 

Power:  Similar  to  the  current  waveforms,  the  power  of  the  ship  and  PEN  are  unchanged.  The 
power  of  each  HESM,  however,  is  one-half  the  power  requirement  of  the  centralized  HESM 
unit.  The  half-power  requirement  of  each  HESM  unit,  however,  requires  that  both  HESM 
units  operate  in  unison  to  support  PEN  charging. 

Energy:  The  fourth  curve  shows  that  the  energy  requirement  of  each  HESM  unit  is 
approximately  one-half  of  what  it  was  for  the  centralized  HESM  unit.  This  results  from 
assuming  one-half  of  the  moment  of  inertia  for  each  distributed  HESM  unit.  The  flywheel’s 
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speed  and  drifts,  however,  are  similar  to  what  was  observed  before.  A  close-up  of  the 
flywheel  speed  and  energy  is  shown  in  Fig.  12. 

The  following  explanations  compare  result  sets  (c)  and  (g),  or  the  right  side  of  Fig.  7 
(centralized)  to  the  right  side  in  Fig.  11  (distributed). 

Voltages:  Similar  to  the  ramped  charging  profile,  the  voltage  of  each  HESM  unit  is  much 
closer  to  the  ring  bus  voltage  than  it  was  before. 

Current:  The  ship  current  remains  unchanged,  which  means  the  combined  current 
contribution  of  both  generators  is  the  same.  This  combined  ship  current,  as  opposed  to  the 
centralized  topology,  splits  into  each  feeder  serving  the  HESM  units  and  the  feeder  serving 
the  PEN.  The  PEN  current,  however,  is  unchanged  as  in  both  centralized  and  distributed 
cases  the  PEN  remains  as  a  single  unit.  Both  the  charging  and  discharging  currents  of  each 
HESM  unit  also  reduced  to  one-half  of  what  was  required  for  the  centralized  HESM  unit. 
Power:  The  power  contributed  by  the  generators  {Pship)  and  the  power  demanded  by  the  PEN 
(Pp/„)  do  not  change.  What  changes  is  the  charging  and  discharging  power  of  the  HESM 
unit.  The  charging  and  discharging  power  of  each  HESM  unit  are  one-half  of  what  they  were 
before. 

Energy:  The  energy  requirement  for  each  HESM  unit  is  approximately  one-half  of  what  it 
was  for  the  centralized  HEMS  unit.  This  result  stems  from  assuming  one-half  the  moment  of 
inertia  for  each  distributed  unit  as  compared  to  the  centralized  HESM  unit.  The  flywheel’s 
speed  and  drifts,  however,  do  not  change.  A  close-up  of  the  flywheel  speed  and  energy  is 
shown  in  Fig.  12. 

7.2.1.2  Comparison  of  ramped  and  hybrid  charging  profdes 

The  following  explanations  compare  the  left  and  right  sides  of  Fig.  11  (result  sets  (e)  and 
(g)). 

Voltages:  The  voltages  are  similar  in  both  cases.  Due  to  the  high  power  peaks  in  Ppfa  during 
the  ramped  charging  profile,  the  voltage  at  the  PEN  terminals  is  slightly  reduced.  This 
reduction,  however,  is  contingent  on  proper  conductor  selection  for  the  PEN  cable  (see  Table 
1). 

Current:  The  ship  current  Chip  and  PEN  current  /p/„  are  similar  in  both  cases.  What  changed 
is  the  peak  charging  and  discharging  HESM  currents.  Referring  to  Ihesmi  (same  as  hesmi),  the 
peak  discharging  current  is  -550  A  or  28  %  of  the  peak  value  of  1 .9  kA  for  the  ramped 
charging.  The  HESM  charging  current  is  the  same  in  both  charging  variants  at  1.35  kA. 

Power:  The  third  curve  set  shows  approximate  power  demands  for  the  ship  in  both  ramped 
and  hybrid  charging  profiles.  The  peak  PEN  charging  power,  in  a  manner  similar  to  the 
centralized  case  (Eig.  7  (a)  and  (c)),  is  also  reduced  from  Ppfr  ~  40  MW  to  Rp/„~23  MW  with 
the  hybrid  charging  profile.  The  hybrid  charging  profile  also  reduces  the  peak  HESM  power 
from  1 1.5  MW  to  8.5  MW  (per  unit). 
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Energy.  The  energy  stored  in  the  HESM  is  shown  in  terms  of  flywheel  speed  (RPM)  on  the 
right-side  axis.  Similar  to  the  centralized  HESM  topology,  flywheel  speed  in  the  hybrid 
charging  profile  reduces  the  ending  swings  by  ~25  %  (100x(631-469)/63 1=25.67). 

7.2.2  Simulation  with  Actual  Sources 

The  results  of  the  simulations  with  actual  sources  are  summarized  in  Fig.  13  (left)  for 
constant  current  charging  and  right  for  hybrid  charging  result  sets  (f)  and  (h).  As  can  be  seen, 
they  reproduce  very  closely  the  results  obtained  with  ideal  sources  (Fig.  11)  so  that  the 
observations  made  above  for  the  case  of  ideal  sources  hold  for  the  case  of  actual  sources  as  well. 
Likewise,  the  comments  already  made  for  the  case  of  a  single  HESM  (Fig.  9)  carry  over  to  the 
present  case  of  two  HESM  (Fig.  11).  The  traces  for  the  physical  quantities  related  to  each  HESM 
are  close  and  overlap. 
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Fig.  11.  Ideal  simulation  results  for  ramped  and  hybrid  PFN  charging  variants  (distributed  HESM).  (Actual  simulation  results  are  shown  in  Fig.  13.) 
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Fig.  12.  Ideal  flywheel  speed  and  energy  for  ramped  and  hybrid  PFN  charging  variants.  (Close-up  of  the  lower  row  of  traces  in  Fig.  11.) 
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Ramped  PFN  Charging  (variant  1) 

vAc^ya/ Simulation  Results 


(h) 

•  Distributed  HESM 

•  Hybrid  PFN  Charging  (variant  2) 

•  vAcfua/Simulation  Results 


Fig.  13.  Actual  simulation  results  for  ramped  and  hybrid  PFN  charging  variants  (distributed  HESM).  (Ideal  simulation  results  were  shown  in  Fig.  11.) 
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8  Simulation  Models 

High  fidelity  simulation  models  were  developed  to  represent  the  two  architectures  and 
evaluate  the  performance  of  the  centralized  and  distributed  energy  storage  modules.  Figure  14 
shows  the  top  level  of  one  of  the  models;  each  major  block  has  operator  inputs  to  characterize 
the  systems  under  simulation. 

The  model  shown  was  developed  with  the  assumption  of  two  ship’s  generators,  one  PFN 
with  pulsed  load,  and  four  HESMs  to  support  the  PFN.  The  modularity  of  the  model  allows  the 
user  to  reposition  all  elements,  generators,  PFN,  and  HESMs,  as  desired  along  the  ship’s  power 
system.  For  example,  the  HESMs  are  shown  in  Figure  14  distributed  around  the  ship’s  power  bus 
but  could  be  all  located  and  connected  at  the  same  point. 

Furthermore,  additional  elements  could  be  added  or  some  of  the  existing  components  could 
be  removed  as  desired  provided  some  basic  control  parameters  are  properly  adjusted  to  account 
for  the  changes.  This  can  easily  be  done  by  following  the  instructions  in  the  data  initialization 
input  file,  also  provided  with  the  model.  For  example  and  for  the  convenience  of  the  Navy,  a 
second  version  of  the  model  shown  in  Figure  14  was  also  provided  with  two  generators,  one 
PFN,  and  two  HESMs.  Thus,  many  reference  models  can  easily  be  generated  from  the  two 
provided  to  the  Navy.  Models  with  up  to  eight  HESMs  have  been  run  successfully,  although 
execution  time  is  strongly  dependent  on  the  complexity  of  the  system  modeled. 

Interesting  issues  that  can  be  analyzed  with  the  models  developed  do  not  only  concern  the 
load  side  of  the  system  but  also  the  source  side.  Thus,  for  example,  one  can  study  the  effect  of 
the  transient  load  at  the  terminals  of  the  ship’s  generators  and  the  resulting  power  quality  on  the 
ship’s  bus  to  determine,  for  example,  if  the  limits  set  forth  in  MIL  STD  1399  are  respected.. 

Likewise,  one  can  monitor  how  the  load  is  distributed  among  the  generators  present  when  a 
transient  occurs.  In  this  latter  case,  it  was  interesting  to  note  that,  depending  on  where  the 
generators,  the  PFN,  and  the  HESMs  were  located  relative  to  each  other,  there  was  a  greater  or 
lesser  likelihood  that  the  generators  would  be  sharing  the  load  somewhat  equally.  In  fact,  there 
were  cases  in  which  a  generator  would  end  up  with  a  disproportionate  fraction  of  the  load 
leading  at  times  even  to  conditions  of  overload.  Since  this  is  a  very  real  concern,  the  models 
were  also  equipped  with  an  optional  power  balancing  circuit  that  tries  to  keep  the  load  equally 
distributed  among  the  existing  generators. 
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Genef  ator  Model 


Slup  Model 


Hybrid  Energy  Storage  Moduli^ 


Pulse  Forming  Network  Module 


Fig.  14.  One  of  the  two  reference  MatLab  Simulink  models  for  the  study  of  the  effect  of  HESMs  and  a  pulsed 

load  on  a  ship’s  power  system. 

Annotated  simulation  models  are  attached  to  this  report  as  Appendix  1 . 
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9  Machine  Characteristics 

The  power  and  energy  delivery  parameters  identified  through  the  simulations  run  sets  ((a) 
through  (h)  in  Fig.  1)  provide  initial  insights  into  the  design  of  the  HESM  flywheel  battery:  the 
combination  of  the  rotating  inertia  and  the  electric  machine  required  to  charge  and  discharge  the 
energy  store.  The  final  design  of  the  rotating  machine  for  the  HESM  will  require  additional 
analysis.  Based  on  the  initial  simulation  results  in  Fig.  7-Fig.  9  (result  sets  (a)  through  (d)),  the 
peak  power  rating  for  the  HESM  ranges  from  approximately  23  MW  for  the  constant 
current/ramped  power  charging  profile  to  about  17  MW  for  the  hybrid  charging  profile.  The 
energy  transfer  energy  from  the  HESM  ranges  from  approximately  28  MJ  (left  side  in  Fig.  8)  for 
the  constant  current/ramped  power  charging  profile  to  about  20  M.T  for  the  hybrid  charging 
profile  (right  side  of  Fig.  8). 

Using  the  ratio  of  peak  power  to  delivered  energy  as  a  figure  of  comparison,  the  constant 
current/ramped  power  profile  has  a  power  to  energy  ratio  of  1 .42  and  the  hybrid  charging  profile 
has  a  power  to  energy  ratio  of  0.85.  In  contrast,  a  flywheel  battery  for  a  typical  vehicular  or  UPS 
application  will  have  power  to  energy  ratios  at  least  an  order  of  magnitude  lower;  for  example, 
two  hybrid-electric  vehicle  flywheel  programs  at  UT  have  power  to  energy  ratios  of  0.02  and 
0.004.  These  figures  indicate  that  the  HESM  design  is  heavily  driven  by  the  power  requirement 
of  the  system  as  opposed  to  the  amount  of  stored  and  delivered  energy.  Electrochemical 
batteries  are  not  well  suited  for  this  type  of  power  and  energy  profile  —  they  are  generally  more 
appropriate  in  applications  with  a  lower  power  to  energy  ratio.  (This  characteristic  may  change 
if  the  HESM  is  used  to  provide  energy  storage  for  other  shipboard  applications  such  as  UPS  for 
single  generator  set  operations.) 

The  HESM  power  requirements  primarily  impact  the  design  of  the  motor/generator  and  the 
power  converter;  however,  the  optimum  HESM  design  must  consider  a  system-level 
optimization  of  all  three  major  elements  of  the  HESM  system  including  the  rotating  inertia,  the 
motor/generator,  and  the  bi-directional  power  converter.  The  following  sections  present  the 
critical  design  issues  for  each  element. 

9.1  Power  Conversion  Equipment 

The  motor/generator  maximum  operating  speed  will  be  limited  by  the  characteristics  of  the 
bi-directional  power  converter  -  the  converter  must  be  able  to  provide  the  required  power  at  the 
fundamental  frequency  to  match  the  speed  and  pole  count  of  the  motor/generator.  Although  it 
increases  the  fundamental  frequency  required  for  a  given  speed,  a  higher  pole  count  is  desirable 
because,  for  a  given  power  rating,  the  power  density  of  the  machine  will  tend  to  increase  as  the 
number  of  poles  increases  [13]. 

In  general,  the  frequency  capability  of  a  power  converter  will  decrease  as  the  power  rating  of 
the  converter  increases.  However,  the  design  can  take  advantage  of  increased  phase  counts  for 
the  driven  rotating  machine,  and  a  higher  power  motor  can  be  driven  using  multiple  smaller 
converters.  For  example,  a  12-phase,  16  MW  motor  can  be  driven  by  an  array  of  four  three- 
phase  converters,  each  with  a  nominal  rating  of  4  MW.  This  approach  also  potentially  offers 
increased  reliability  -  if  one  converter  module  fails,  the  motor  can  continue  operation  at  reduced 
power.  Moreover,  this  approach  also  applies  to  architectures  based  on  dc  or  ac  distribution 


31 


systems  since  the  intermediate  dc  link  of  ac-ac  converters  provides  the  compensating  point  for 
the  phase  shifts  among  the  parallel  three-phase  units. 

Although  silicon-carbide  and  gallium-nitride  semiconductor  technologies  offers  the  promise 
of  significant  increases  in  converter  performance,  these  technologies  are  not  yet  mature  enough 
for  critical  DOD  applications.  For  the  purposes  of  this  study,  it  is  assumed  that  the  design  will  be 
limited  to  a  maximum  fundamental  operating  frequency  of  500  Hz  and  that  multiple  converter 
modules  will  be  used  to  drive  higher  power  machines.  Given  the  anticipated  adoption  of  wide 
band  gap  semiconductor  switching  at  higher  speeds,  this  should  be  considered  a  lower  bound  on 
the  future  capability. 

9.2  Machine  Considerations 

The  inertia  of  the  HESM  flywheel  rotor  will  define  the  maximum  stored  energy  at  a  given 
angular  velocity: 

W  =  ^Ja)^  (8.1) 

where 

W  =  energy  stored  in  flywheel  mass  [Joules] 

J  =  moment  of  inertia  [kg-m^] 
o)  =  flywheel  speed  [rad/s]. 

Assuming  a  depth  of  discharge  of  50%  of  maximum  speed,  the  energy  transferred  from  the 
flywheel  storage  is  equal  to  75%  of  the  maximum  stored  energy.  This  also  means  that,  knowing 
the  minimum  amount  of  energy  that  must  be  transferred,  we  must  design  the  machine  such  that 
the  stored  energy  is  at  least  33%  larger  than  the  delivered  energy.  This  figure  defines  the 
minimum  energy  storage  at  maximum  operating  speed;  the  actual  energy  storage  will  likely  be 
larger  than  this  value.  Depths  of  discharge  below  this  value  offer  diminishing  returns  and 
increase  the  challenges  for  the  motor/generator  to  make  required  voltage  at  the  lower  speeds. 

While  energy  storage  flywheels  allow  independent  selection  of  energy  and  power,  the  design 
of  the  HESM  inertia  cannot  be  defined  without  consideration  for  the  characteristics  of  the 
motor/generator  and  power  converter.  The  squared  relationship  of  energy  with  angular  velocity 
drives  the  design  towards  the  highest  practical  rotational  speed  to  maximize  energy  density; 
however,  the  relatively  high  peak  power  of  1 7  to  40  MW  effectively  limits  the  minimum 
physical  size  and  maximum  operating  speed  of  the  motor/generator.  (The  peak  power  values  are 
of  relatively  short  duration  so  the  designs  can  be  more  aggressive  than  machines  rated  for 
continuous  duty  at  these  power  levels.)  This  can  be  illustrated  by  consideration  of  the  maximum 
practical  tip  speed  of  the  flywheel  and  motor/generator  and  the  effective  torque  generated  per 
unit  surface  area  or  unit  volume  of  the  motor/generator. 

Tip  speed  simply  refers  to  the  product  of  angular  velocity  and  outer  radius  of  the  rotor.  UT 
has  demonstrated  high  performance  graphite  epoxy  composite  flywheels  at  tips  speeds  in  excess 
of  1,350  m/s  and  small  devices  have  reached  tip  speeds  of  over  2,000  m/s.'*  In  contrast,  the 
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weight  of  permanent  magnets  or  windings  limits  the  maximum  practical  tip  speed  for  an  electric 
machine  to  approximately  3-400  ni/s.  For  flywheel  applications,  the  speed  decreases  as  you 
extract  energy  from  the  rotor  so  motor/generator  must  be  capable  of  the  required  performance  at 
the  minimum  speed  of  the  machine.  These  constraints,  in  conjunction  with  electromagnetic  and 
rotor  dynamic  considerations,  define  the  size  of  the  motor/generator. 

Torque  generation  in  an  electric  machine  is  described  by  the  following  relationship: 

TccBxA^xLxD^  (8.2) 

where 

r=  torque  [N-m] 

B  =  air  gap  flux  density  [Wb/m  ] 

As  =  stator  line  current  density  [A/m] 

L  =  active  length  [m] 

D  =  air  gap  diameter  [m]. 

The  air  gap  flux  density  is  set  by  the  B-H  (magnetic  flux  density  vs.  intensity)  characteristics 
of  the  iron  in  the  rotor  and  stator  and  by  the  available  magneto-motive  force  (MMF)  of  the 
excitation  source  (the  strength  of  permanent  magnets  or  the  allowable  current  or  current  density 
in  the  field  windings).  In  slotted  structures,  the  flux  density  in  the  teeth  is  approximately  twice 
that  of  the  air  gap  so  for  a  typical  saturation  limit  of  1 .8  T,  the  peak  field  in  the  air  gap  is  limited 
to  approximately  0.9  Tpeak,  or  about  0.57  T  rms. 

The  diameter  of  the  air  gap  is  maximized  to  take  advantage  of  the  squared  relationship  but  is 
limited  by  the  practical  consideration  of  the  3-400  m/s  tip  speed  limit.  The  length  of  the  rotor  is 
effectively  constrained  by  rotor  dynamic  considerations — as  the  rotor  gets  smaller  in  diameter 
and  longer,  its  stiffness  is  not  adequate  to  allow  sub-critical  operation  at  high  speed.  This  is 
particularly  important  for  flywheel  applications  where  the  speed  is  changing  during  every 
charge/discharge  cycle.  Supercritical  rotors  may  require  negotiation  of  the  critical  speeds  on 
every  cycle. 


The  torque  being  generated  at  the  rotor  surface  can  be  linked  to  the  volume  of  the  rotor  [14]: 

(8,3) 

TTr  I 


where 

TRV  =  torque  per  unit  rotor  volume  (kN-m/m^) 
T  =  torque  generated  (kN-m) 

r  —  rotor  radius  (m) 

/  =  active  length  (m). 


If  one  considers  the  tangential  force  generated  by  the  J  x  .S  interaction  of  the  current  and 
magnetic  field  in  the  air  gap  and  the  surface  area  of  the  rotor,  an  effective  electromagnetic  “shear 
stress”  can  be  calculated.  Referring  to  the  “shear  stress”  as  sigma,  a,  this  parameter  can  be 
related  to  the  torque  generated  per  unit  rotor  volume,  TRV  as: 

TRV  =  2a  (8.4) 

where 
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a  =  shear  stress  in  kN/m^. 


High  performance,  liquid  cooled  electric  machines  typically  demonstrate  shear  stresses  on 
the  order  of  70  -100  kPa  and  TRV  values  between  130  and  220  kN-m/m^. 


Using  these  relationships,  one  can  begin  to  estimate  the  size  of  the  motor  generator  as  a 
function  of  minimum  rotor  speed.  Thus,  from  rotor  angular  velocity  and  for  a  given  tip  speed  the 
rotor  radius  can  be  calculated. 

Furthermore,  recalling  that: 

P  =  T -CO  or  T=—  (8.5) 

CO 

where 

P  =  power  [W] 

7"=  torque  [N-m] 

0)  =  angular  velocity  [rad/s]. 


Fig.  15  shows  the  calculated  torque  for  a  variety  of  operating  speeds  for  a  single  HESM  as  well 
as  for  multiple  units. 


HESM  Motor/Generator  Torque  vs.  Speed 


Fig.  15.  Torque-speed  curves  for  motor/generator  sets  of  different  size. 


Using  these  torque  values  and  a  representative  value  for  effective  shear  stress,  one  can 
calculate  the  required  surface  area  of  the  motor/generator  rotor  and  from  this,  since  the  radius  is 
already  known,  the  rotor  length  can  be  derived.  Then  the  rotor  moment  of  inertia  and  the  energy 
stored  in  the  rotor  can  be  calculated. 


Proceeding  in  this  fashion,  one  can  make  the  following  preliminary  estimates: 

1 .  A  single  HESM  unit  with  a  rotor  radius  of  40.9  cm,  a  rotor  length  of  25.6  cm,  spinning  at 
7,000  RPM  (467  Hz  from  the  converter  for  an  8-pole  machine)  can  store  a  maximum  of 
23.8  MJ  and  operates  at  17  MW. 
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2.  Two  HESM  units  with  a  rotor  radius  of  40.9  cm,  a  rotor  length  of  12.8  cm,  spinning  at 
7,000  RPM  (467  Hz  from  the  converter  for  an  8-pole  machine)  can  store  a  maximum  of 

1 1.9  MJ  and  operate  at  8.5  MW.  Thus  the  combined  performance  of  the  two  units  is 
equal  to  that  of  the  single  HESM. 

Another  possibility  is: 

1.  A  single  HESM  unit  with  a  rotor  radius  of  23.9  cm,  a  rotor  length  of  43.9  cm,  spinning  at 
12,000  RPM  (400  Hz  from  the  converter  for  a  4-pole  machine)  can  store  a  maximum  of 

13.9  MJ  and  operates  at  17  MW. 

2.  Two  HESM  units  with  a  rotor  radius  of  23.9  cm,  a  rotor  length  of  22.0  cm,  spinning  at 

1 2,000  RPM  (400  Hz  from  the  converter  for  a  4-pole  machine)  can  store  a  maximum  of 
7.0  MJ  and  operate  at  8.5  MW.  Once  again,  the  combined  performance  of  the  two  units  is 
equal  to  that  of  the  single  HESM. 

Likewise,  more  combinations  can  be  considered.  It  can  be  noticed  that  the  performance  of  the 
first  system  so  estimated  comes  very  close  to  fulfilling  all  requirements  for  our  HESM  system 
operated  in  the  hybrid  charging  mode.  Thus,  essentially  no  external  flywheel  is  necessary  in 
addition  to  the  energy  stored  by  the  HESM  rotor  itself,  provided  the  rotor  is  made  a  little  longer 
than  required  to  deliver  the  required  power.  This  is  quite  feasible,  since  the  aspect  ratio 
(length/diameter)  of  the  rotors  calculated  above  is  substantially  less  than  1.0  (0.313  and  0.157 
respectively). 

Conversely,  the  second  system  estimated  above,  will  need  additional  inertia  to  fulfill  the 
requirements  of  our  HESM.  In  this  case,  however,  the  aspect  ratio  of  the  machine’s  rotors 
calculated  above  is  closer  to  1 .0  (0.92  and  0.46  respectively)  and,  thus,  will  likely  lead  to  a  better 
design  from  the  standpoint  of  volumetric  power  density. 

9.3  Machine  Sizing  Study 

The  goal  of  this  task  is  to  quantify  the  trade-offs  for  the  number  of  flywheels  relative  to  the 
total  size  and  mass  devoted  to  flywheel  energy  storage  using  detailed  sizing  studies  based  on 
previous  UT-CEM  flywheel  programs. 

The  following  Figure  16  summarizes  the  various  flywheel  designs  developed  at  UT-CEM 
and  highlights  the  two  models  used  for  this  initial  sizing  study. 
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D«ta  it  For  omelai  Govemmant  U«a  Only 

Machina  Name 

Peak  Power 
(kW) 

Continuous 
Power  (kW) 

Rotor 

Energy 

(MJ) 

Rotor 

RPM 

Rotor  Design 

Tip  Speed 
at  Fun 
Energy 
(m/sj 

oeptn  of 
Dtsctiarge 
(energy,  %) 

Mass  of 

Powar 

Electronics 

Mass  of 
Flywtieat 
Moduia  (IncI 
motoT'gan] 
(kg) 

Power 
Electronics 
Dimensions  (m) 

Flywheel 

Module 

Dlametar 

(m) 

Rywheel 
Module  Length 
(m) 

Project 

8tatus 

Comments 

EMAUe 

60,000 

n/a 

117.0 

4,200 

Partially  integratad, 
cylindrical,  masa 
loaded,  non*composlte 

240 

75% 

19,273 

1.75 

3.10 

Built 

EMALS  Is  a  ■yncltfonoua  meelWM  (ttt«  sntrgv  storwl  wHh  an  100  kW  aacHar  and  a  1.2  MW 
Inductton  motor  to  laeecatarata  ttw  mechlfM  altar  Via  diacltarga,  all  In  line  on  the  aama  ahaft. 
la  a  pula*  powar  machltta.  Ttw  alactronica  ware  not  part  of  CEM'a  afiort 

ALPS 

4,000 

2,000 

478,8 

15,000 

Non-lntagratad, 
cylindrlcai,  composite 

900 

75% 

3,000 

13,400 

2  X  2.3  X  1.4 

1.84 

1J8 

Built 

ALPS  la  a  non4ntagrat«d  machina,  wtth  tfia  rnotor  ganarator  aaparata  lorm  tha  anargy  atoraga 
rotor. 

Trentft  Bus 

250 

110 

7.2 

40,000 

Partially  Integrated, 
cyilndilcai,  composite 

935 

75% 

250 

0.7x0.5x0.2 

0.70 

0.75 

Bum 

MaeMtia  waa  capaWa  ol  2M  kW.  but  powar  elactrorrlea  wara  only  eapabla  oMM  kW.  MacNna 
watgirt  doaa  not  mduda  contalnmant  or  gicnbaia. 

CMPS 

10,000 

350 

25.2 

20,000 

Fully  integrated  inside 
out  mass  ioaded, 
composite 

800 

75% 

112 

480 

Designed  ; 
ttifough 
COR 

Duel  purpoaa  machina.  wHh  5UW  rebrtg  ter  pulaa  waapona  and  560  kW  tor  hybrid  atacOic 
combat  yaMM  toad  levaSng. 

FCSS 

5.0 

3.8 

13.0 

53,000 

PardaRy  imegrated, 
cylindrical,  compoaRa 

920 

88% 

115 

Designed 
thru  COR 

Daslgn  tecus  on  high  offictancy  (>S4%  turn  around  alSclancy  ter  Sywhaat] 

Small  Cal 

864,000 

n/a 

g.o 

25,000 

ComposKe,  cylindrical, 
mesa  loaded 

502 

850 

0.40 

0.43 

Built 

Pulaa  powar  machine 

Cannon  Cal 

2.400.000 

n/a 

30.8 

12,000 

Composite,  cyilndrlcal, 
mass  loadad 

458 

50% 

865 

0.73 

1.00 

Built 

Pulaa  powar  meehina 

Sub  Scale 

3,000,000 

n/a 

23.0 

12,000 

Composlta  arbor, 
composite  rotor 

380 

2,400 

0.80 

1.00 

Bum 

Pulaa  powar  machina 

AIR 

73,000 

n/a 

6,2 

8,400 

AB  Iron  Rotor 

Built 

Pulaa  powar  machina 

iron  Cora 

1,000,000 

n/a 

37.8 

4,700 

iron 

200 

10,000 

0J10 

1.00 

Built 

Pulaa  powar  tnaclXna 

HVHP 

125,000 

n/a 

3.25 

8.800 

Aluminum 

220 

0.84 

Built 

Pulaa  powar  machina:  homopolar  machine 

BHPC 

00,000 

n/a 

10.0 

8,073 

Iron 

181 

8,000 

0.57 

0.59 

Built 

Pulse  powar  machina:  homopolar  machina 

Fig.  16  UT-CEM  flywheels  and  the  two  models  selected  for  the  initial  sizing  study.  The  BUS  flywheel  is  a  partially  integrated  design  and  the  CHPS 

machine  is  a  fully  integrated  design. 
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Fig,  17  BUS  and  CHPS  flywheel  batteries  used  in  initial  sizing  study. 

Matlab  Simulink  models  were  developed  as  a  design  tool  to  calculate  the  size  and  weight  of 
flywheel  rotors  over  a  range  of  energy  and  power  ratings  for  the  two  topologies  presented  above. 
The  rotor  sizing  tool  was  validated  against  the  actual  rotor  designs  for  the  two  programs.  Figures 
18  and  19  show  representative  models  of  the  two  flywheels  used  for  the  validations. 


Upper  Back  Up  Bearing  Interface 
Radial  Mag  Bearing  Interface 
Upper  Hub/Shaft 
Permanent  Magnet  (PM)  Cartridge 
Lower  Hub/Shaft 
Composite  Energy  Storage  Rings 


Lower  Back  Up  Bearing  Interface 


Fig.  18  Rotor  sizing  tool  for  partially  integrated  flywheel  design. 
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Based  on  the  results  of  the  initial  evaluations  using  the  rotor  design  tool,  a  more  detailed 
study  was  conducted  using  the  CHPS  baseline  design  topology.  There  are  two  strong  indications 
that  the  CHPS  flywheel  design  can  provide  an  approach  to  meeting  the  baseline  HESM 
requirements  with  only  two  machines  that  each  can  fit  through  a  26  inch  (0.66  m)  hatch. 

•  Use  TRV  scaling  approach  to  show  that  CHPS  topology  and  envelop  is  within  the  state- 
of-the-art  (referred  to  as  ''TRV scaling'  approach). 

•  Use  CHPS  program  plots  to  provide  a  direct  path  to  meet  HESM  requirements  by  adding 
an  extra  energy  storage  ring,  which  increases  flywheel  housing  OD  from  22  inches  to  24 
inches  (referred  to  as  "added  composite  rin^'  approach. 

For  comparison  with  the  HESM  requirements  based  on  hybrid  charging  profile: 

•  HESM  required  power  (occurs  during  flywheel  charging):  17  MW  (8.5MW  per 
machine) 

•  HESM  rms  power  required:  9.3  (4.65  MW_rms  per  machine) 

•  HESM  required  delivered  energy:  20MJ  (lOMJ  per  machine) 

Baseline  CHPS  design  is  good  match  with  ~25%  depth  of  discharge: 

•  Wmax=20,000RPM;wn,in  =  14,859RPM 

•  Power  at  Wmin  =  8.52MW 

•  E  wmax  =  24.59  MJ;  E  wmin  =  13.57  MJ;  Energy  Delivered  =11  MJ 

The  CHPS  machine  was  originally  designed  for  operation  with  a  resonant  converter  that 
allowed  for  higher  performance  than  is  projected  by  the  TRV  evaluation.  The  system  is  very 
attractive  for  mobile  or  space  constrained  applications  due  to  the  high  volumetric  energy  and 
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power  density  of  the  custom  power  supply.  Figure  20  shows  the  components  of  the  CHPS  power 
supply  to  scale;  the  OD  of  the  CHPS  flywheel  battery  will  fit  through  a  26”  hatchway. 


Fig.  20  CHPS  power  supply  components  to  scale. 


Output  power  is  not  linear  with  rotor  speed  due  to  resonant  output  circuit  so,  for  the  second 
design  evaluation,  the  HESM  was  projected  using  data  from  the  power/speed  curves  including  a 
custom  resonant  converter.  Figure  21  shows  the  power  and  efficiency  performance  as  a  function 
of  stored  energy  (rpm),  into  a  resistive  load  including  saturation  effects. 


25  24  23  22  21  20  19  18  17  16  15  14 

Stored  Energy  (MW) 

Fig.  21  Power  and  efficiency  as  a  function  of  stored  energy  into  a  resistive  load  including  saturation  effects. 

For  this  design,  an  additional  flywheel  ring  was  also  incorporated  onto  the  rotor  to  increase 
the  inertia  and  stored  energy.  Due  to  the  high  tip  speed  of  this  application,  the  composite  ring  is 
constructed  of  high  modulus  TIOOOG  graphite/epoxy  composite.  The  strains  in  the  outer 
composite  ring  are  well  within  demonstrated  performance  and  have  significant  margin  at  the 
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operating  point  required  for  the  HESM  application.  This  implies  that  additional  gains  in  stored 
energy  are  possible,  albeit  at  the  expense  of  reduced  structural  margin  and  increased  risk.  Table 
2  summarizes  the  strain  and  stress  characteristics  of  the  added  TIOOOG  flywheel  ring. 


Table  2.  Stress  and  strain  of  TIOOOG  flywheel  ring  over  operating  temperature  range. 


75°F 

0°F 

225°F 

Factor  of ' 
Safety 

0  rpm 

1 5  krpm 

20  krpm 

0  rpm 

15  krpm 

20  krpm 

0  rpm 

15  krpm 

20  krpm 

FI  radial  compression 

-13.3 

-14.6 
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-11.1 

-12.4 

-16.7 

-17.9 
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2.2 

axial  tension 

3.4 

3.4 

3.3 

3.6 

3.5 

3.4 

2,2 

2.1 

2 

3.9 

rz  shear 

2.6 

2.8 

3 

2.1 

2.3 

2.5 

3.9 

4,1  i 

t  4.3 

1.9 

hoop  strain 

0.0019 

0.0037 

0.0061 

0.0018 

0.0037 

0.0061 

0.0024 

0.0041 

‘  0.0061 

2.3 

F2  radial  compression 

-5.7 

-5.5 

-6.6 

-4.7 

-4.6 

-5.8 

-7.7 

-7.5 

-7.8 

5.8 

axial  tension 
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2.2 

2.7 

2.8 

2.9 

3.4 

2.4 

2 

1.7 

4.1 

rz  shear 

1.8 

2 

2.1 

1.4 

1.6 

1.9 

2.5 

2.7  1 

1  2.8 

2.9 

hoop  strain 

0.0021 

0.0035 

0.0052 

0.0019 

0.0033 

0.005 

0.0027 

0.0041 

'  0.0056 

2,5 

F3  radial  compression 

-0.81 

-0.72 

-0.78 

-0.71 

-0.62 

-0.68 

-1.1 

-1 

-0.98 

45.9 

axial  tension 

3.8 

3.3 

3.7 

3.8 

3.7 

4.3 

3.8 

3.3 

2.9 

3.3 

rz  shear 

0.32 

0.36 

0.38 

0.26 

0.29 

0.36 

0.4 

0.5 

0.5 

16.0 

hoop  strain 

0.0024 

0.0036 

0.0051 

0.0021 

0.0033 

0.0047 

0.0032 

0.0044 

0.0057 

2-5  ! 

Radial  compression 

45 

ksi 

Note:  All  strength  and  stress  units  are  KSI 

Axial  tension 

14 

ksi 

rz  shear 

8 

ksi* 

Hoop  strain 

1.40% 

ksi 

(from  elevated  temperature  tests  at  275  ®F) 

*  from  earlier  database,  needs  to  be  verified 


The  CHPS  flywheel  topology  provides  an  effective,  high  power  density  energy  storage 
module  capable  of  addressing  the  performance  requirements  of  the  PEN  charging  application. 
Tables  3  and  4  show  the  rotor  design  and  machine  envelope  dimensions  generated  using  the  rotor 
design  tool. 
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Table  3.  Baseline  CHPS  sizing  using  TRV  calculations. 


MATLAB  RESULTS 


Rotor 


B/U  Brg  1x2} 

Combo  MB 

Combo  MB  KSG 

Radial  MB 

Radial  MB  KS6 

Composite  Rings 

PM  Back  Iron* 

PM  Cartridge 
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Low  Back  Iron* 

Top  Back  Iron* 
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Mass  (kg) 

14.2 

14.S 
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17.0 

3.3 
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49.S 

36.7 

2.0 
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4.2SE'01 

0.31 
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0.11 

6.S6 

1.61 

1.00 
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Jt  Ikg-m*!) 

2.46E4O0 

1.2S 

1,16 
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11.02 

1.22 

0.81 
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0.27 

0.55 

21.19 

Stored  Energy  (MJ) 

0.93 

0.67 

0.8S 

1.06 

0.24 

14.39 

3.54 

119 

0.11 

0.19 

0.42 

24.59 

‘Note:  Rotating  Back  Iron  is  divided  into  PM  Back  Iron,  Lower  Back  Iror>  and  Upper  Back  iron  sections  for  computational  puiposes. 


Stator  ar>d  Miscellar>eous  Performance  Results 


Vac  Hsg  Mass 

Vac  Hsg  OD  1  VacHsgOAL  I  Shaft  Mass  |  MG  Mass  1  Total  BrgMass  '  Stator  Total  Mass  |  SbearUmPwrMaKfM])  I  ShearUmPwrMinlMJ)  1  Rotor  CM 

Jp/Jt  1  MinRPM  1  Discharge  Depth 

Del  Energy  IMJJ 

97.2 

0,550  1  0.986  1  84.1  511  |  71.2  304.6  11.5  8.5  |  0.427 

0.53  1  14,859  26% 

11.02 

Housing  OD  (in)  =  21.7  Mass  (kg)  =  59B  Volume  (liter)  =  234 


Table  4.  Optimized  CHPS  design  with  additional  T1 OOOG  flywheel  ring. 


MAT 

U5  RESULTS 

Rotor 

B/UBrg(x2) 

Combo  MB 

Combo  MB  h5G 

Radial  MB 

Radial  MB  HSG 

Composite  Rings 

PM  Back  Iron* 

PM  Cartridge 

PM  Sleeve 

Low  Back  iron* 

Top  Back  Iron* 
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Ma«(k|) 

14.2 

14.5 

12.9 
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33 
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36.7 

2.0 

2.5 
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0.31 

0.39 

0.48 
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10.96 

1.61 
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0.05 
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15.60 
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146£tOO 

U6 
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16.S6 
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O.Sl 

0.04 
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Stored  Energy  (MJ] 
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0.67 

025 

1.06 

0.24 
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334 
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0.11 

0.19 

0.42 
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•Note:  Rota 

itirg  Back  Iron  is  divided  into  PM  Back  Ir 

on.  Lower  Back  Iron  and  Upper  Sack  Iron  actions  for  compL'tationa!  purposes. 

Stator  and  MisceBaneotis  Performance  Results 

Vac  K 94  Mass  VachhOD 

VacHsi  OAL 

5h8ft  Mass 

MG  Mass 

Total  Brg  Mass  I  Statoc  Total  Mass 

SheaTihnPwrMButMJ]  1  ShearliniPwTMinlMJI  I  Rotor  CM 

Jp/Jt  Min  RPM 

Dischafbe  Oejth  Del  bfergy  (MJ) 

llljO  0.601 

09S6  ^ 

£4.1 
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71.2  1  318.4 

115  1  9.6  1  0426 

0.5S  !  16730 

m  1  1078 

Housing  OD  (in)  =  23.7  Mass  (kg)  =  669  Volume  (liter)  =  280 
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10  ALPS  and  MPM  Program  Review 

The  final  task  was  to  assess  the  potential  for  crossover  technology  from  two  earlier 
University  of  Texas  programs:  the  Advanced  Locomotive  Propulsion  System  (ALPS)  Program 
and  the  Megawatt  Power  Module  for  Ship  Service  Program.  These  both  required  modeling  and 
simulation  and  the  development  of  control  algorithms  for  management  of  stored  energy  in 
isolated  power  systems  using  small  gas  turbine  (~3  MW)  prime  movers.  ALPS  Program 


The  goal  of  the  ALPS  program  was  to  demonstrate  the  feasibility  of  a  hybrid  electric 
propulsion  system  for  a  self-powered  locomotive  for  150  mph  high  speed  passenger  rail.  Prime 
power  for  the  locomotive  was  provided  by  a  3  MW  gas  turbine  directly  driving  a  high  speed 
alternator.  A  480  MJ  energy  storage  flywheel  was  used  to  load  level  the  gas  turbine,  recover 
braking  energy,  and  provide  supplemental  power  for  acceleration  and  negotiation  of  grades. 
Figure  22  shows  a  simplified  block  diagram  of  the  ALPS  system;  highlighted  elements  were 
developed  and  demonstrated  during  the  program. 
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Fig.  22  ALPS  system  block  diagram. 


Two  gas  turbine  prime  movers  were  studied  during  the  program,  the  Honeywell  TF40  and 
the  Pratt  &  Whitney  ST-150.  Both  engines  produce  3  MW  at  ISO  conditions  with  a  power 
turbine  speed  of  15,000  rpm.  The  engine  was  directly  coupled  to  an  8-pole  wound  field 
synchronous  generator  developed  in  conjunction  with  AlliedSignal/  Honeywell.  Output  of  the 
high  speed  alternator  was  converted  to  dc  by  a  passive  diode  rectifier  for  powering  the 
locomotives  ~2,000  Vdc  bus.  The  system  also  featured  a  2  MW,  480  MJ  energy  storage  flywheel 
which  operated  between  7,500  and  15,000  rpm.  The  flywheel  included  a  2  MW  induction 
motor/generator  to  provide  the  electro-mechanical  energy  conversion.  The  flywheel  was 
interfaced  to  the  dc  bus  by  a  custom  soft-switching  power  converter  of  the  Auxiliary  Resonant 
Commutated  Pole  (ARCP)  topology. 
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A  top  level  block  diagram  of  the  system  controls  is  shown  in  Figure  23.  The  ALPS  system 
was  designed  for  integration  with  a  locomotive  supervisory  controller  and  much  of  the  overall 
energy  management  functionality  was  outside  of  the  UT-CEM  scope.  The  interface  controller 
provided  a  link  between  the  locomotive  supervisory  controller  and  the  flywheel  and  gas  turbine 
controllers.  Unlike  most  power  system  applications,  because  train  and  detailed  route  information 
are  known  the  demands  on  the  power  and  energy  storage  systems  can  be  defined  a  priori.  This 
situation  is  further  complicated  by  the  requirement  to  replicate  the  existing  8-notch  throttle 
control  architecture  typical  to  locomotives.  This  hybrid  power  system  was  successfully 
demonstrated  in  the  UT-CEM  laboratory. 


10.1  Megawatt  Power  Module  for  Ship  Service 

The  goal  of  the  Megawatt  Power  Module  (MPM)  for  Ship  Service  program  was  to  explore 
the  use  of  high  speed  generation  and  flywheel  energy  storage  to  improve  the  efficiency  of  the 
DDG5 1  ship  service  electric  power  distribution  system  by  enabling  operation  with  a  single 
generator  in  non-threat  environments.  UT-CEM  partnered  with  Rolls-Royce  for  this  program 

To  enable  backfit  onto  existing  platforms,  the  new  system  fits  in  the  same  volume  and 
footprint  as  the  existing  AG9140  turbine  generator  sets.  Concept  of  operations  (CONOPS)  for 
the  DDG5 1  calls  for  two  ship  service  generators  to  be  online  at  all  times  to  prevent  loss  of  power 
(dark  ship)  in  the  event  of  a  failure  of  a  single  unit.  Although  it  contributes  to  higher  reliability, 
the  poor  part-load  efficiency  of  gas  turbines  makes  this  mode  of  operation  highly  inefficient. 
Using  direct  drive  high  speed  generators  both  increased  power  density  and  eliminated  the 
gearbox  which  freed  sufficient  volume  to  incorporate  multiple  energy  storage  flywheels.  The 
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energy  store  provides  un- interruptible  power  (UPS)  for  a  duration  sufficient  for  multiple  start 
attempts  on  a  second  “stand-by”  turbine  generator  set.  The  system  provides  2.5  MW  for  ten 
minutes.  Operation  on  a  single  generator  moved  the  typical  DDG5 1  electric  load  closer  to  the 
engine’s  rated  load,  greatly  improving  specific  fuel  consumption.  The  improved  performance  of 
the  twin  shaft  engine  and  elimination  of  the  gearbox  also  contributed  to  efficiency.  Figure  24  is  a 
rendering  of  one  variant  explored  during  the  program. 
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Fig.  24  General  arrangement  of  one  variant  of  the  Megawatt  Power  Module. 


Figure  25  is  a  block  diagram  of  the  MPM  system  showing  the  baseline  architecture;  several 
variants  were  explored  with  a  range  of  flywheel  and  motor/generator  ratings. 
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Fig.  25  Block  diagram  of  the  MPM  system  —  two  flywheels  per  skid  configuration. 


More  extensive  modeling  and  simulation  was  conducted  for  the  MPM  program  including  the 
development  of  preliminary  control  algorithms  to  manage  the  energy  storage  functions.  Figure 
26  is  an  overview  of  the  MatLab  Simulink  model  of  the  MPM  system  for  a  two  flywheel 
configuration  showing  some  of  the  major  block  functionality. 
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Active  FW 


Of  particular  interest  for  the  HESM  program  is  the  control  algorithm  used  to  manage  the 
stored  energy.  Figure  27  is  a  top  level  flowchart  of  the  control  algorithm  used  in  the  MPM 
modeling. 
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Inputs  from  Sensors  Distributed  in  System 
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Fig.  27  Flowchart  of  the  energy  storage  control  methodology  for  the  MPM  program  modeling. 


11  Summary  and  Conclusions 

The  design  and  integration  of  a  Hybrid  Energy  Storage  Module  (HESM)  to  mitigate  the 
impact  of  transient  (pulsed)  loads  on  the  ship’s  distribution  power  system  has  been  studied.  The 
pulsed  load  was  assumed  to  be  a  high-energy  capacitor-based  Pulse  Forming  Network  (PFN) 
operated  with  a  repetition  rate  of  12  charge/discharge  cycles  per  minute. 

Two  topologies  for  the  HESM  placement  within  the  ship’s  power  system  were  considered: 
centralized  (single  HESM  and  PFN)  and  distributed  (multiple  HESMs  and  PFNs).  Furthermore, 
each  topology  was  evaluated  using  two  PFN  charging  profiles:  constant  current  charging  and 
hybrid  charging.  Each  PFN  charging  profile  was  then  simulated  first  using  controlled  ideal 
sources  for  the  HESM  and  the  PFN,  and  then  using  actual  models  for  the  power  electronic 
converters  generators  and  motors  used  within  the  HESM  and  the  ship’s  power  system. 

The  control  of  the  FIESM  was  accomplished  by  means  of  a  derived  voltage  reference  signal 
(6.4)  that  served  two  purposes:  to  control  the  power  flow  in  and  out  of  the  HESM  controller,  and 
to  maintain  constant  ship  current.  The  simulation  results  validated  this  equation  as  a  proper  way 
to  control  the  HESM’s  bi-directional  converter. 

Also,  an  amp-hour  analysis  suggested  approximate  values  for  Thip  in  each  charging  variant, 
where  this  value  remains  constant  during  PFN  and  HESM  charging  and  discharging.  An 
improper  selection  for  this  value  would  result  in  flywheel  speed  drifts.  These  drifts  were  indeed 
observed  in  tbe  flywheel  speed  simulation  results  and  stem  from  neglecting  the  dc  voltage  ripple 
and  unpredictability  of  the  exact  PFN  peak  charging  current — the  latter  being  a  function  of  its 
terminal  voltage.  The  drift  can  be  eliminated  by  tracking  flywheel  speed  and  adjusting  the  value 
of  Iship  in  (6.4)  or  by  using  the  frequency  control  of  the  power  converter. 

The  charge  and  discharge  power  levels  of  the  HESM  unit(s)  affect  the  size  of  the  bi¬ 
directional  converter.  Although  there  are  many  common  power  conversion  elements,  the  weight 
and  volume  of  the  power  conversion  equipment  is  potentially  a  strong  discriminator  between  the 
various  architectures  and  needs  to  be  addressed  consistently  to  ensure  accurate  comparisons.  The 
ESRDC  developed  basic  volumetric  and  gravimetric  power  densities  for  different  types  of  power 
electronic  conversion  technologies  (e.g.  rectifier,  inverter,  dc-dc  converter)  to  allow  for  a 
common  and  consistent  sizing  approach  for  the  various  power  converters  found  in  the  different 
architectures  presented  herein.  This  approach  should  enable  comparisons  between  the  weight  and 
volume  of  the  power  conversion  equipment  for  different  architectures.  At  this  stage  of  the  study, 
calculation  of  the  weight  and  volume  of  the  equipment  is  based  on  the  power  requirements 
without  attempting  to  identify  commercially  available  devices  or  scale  up  to  the  next  typical 
rating. 

The  work  presented  in  this  report  suggests  the  exclusion  of  the  constant  current  charging 
profiles  since  its  higher  peak  PFN  charging  power  requirements  impose  stricter  specifications  on 
the  PFN  charging  power  supplies  and  HESM  unit  which  reflect  on  the  construction  size  and  cost. 
The  work  showed  that  an  architecture  based  on  a  HESM  concept  to  minimize  the  impact  of 
pulsed  power  loads  on  a  ship’s  power  system  is  feasible  and  very  effective.  Furthermore  it 
confirms  that  the  HESM  can  be  used  as  an  additional  support  element  for  the  ship’s  power 
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system  when  not  needed  for  pulsed  load  charging  (e.g.  UPS  duty),  thus  making  it  a  very  versatile 
component  for  a  variety  of  tasks. 
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12  APPENDIX  1.  HESM  Simulink  Models 


The  following  three  items  listed  in  this  Appendix  have  been  delivered  to  Tracy  Hannon  of 
NAVSEAon  May  1,2013. 

Input  Data  File  HESMm_PFNl_Rl 


U  HESMm_PFNl_Rl.m 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 
11 
12 

13 

14 

15 

16  - 


%  SHIP  WITH  HYBRID  ENERGY  STORAGE  MODULES 
%  Created  by  The  Center  for  Electromechanics  of  The  University  of  Texas  for 
%  The  Office  of  Naval  Research  and  The  Naval  Surface  Warfare  Center 
%  April  30,  2013 
% 

%  Input  items  needed  for  the  simulation  of  a  Ship  Power  System  with  a 
%  number  m  of  Hybrid  Energy  Storage  Modules  and  a  single  Pulse  Forming 
%  Netv/ork, 

% 

make  sure  to  enter  all  required  parameters,  either  directly  or  via 
input  menus,  in  all  blocks  or  elements  with  colored  background 
(not  white)  at  ALL  levels  (remerriber  to  "Look  Under  Mask”)  . 

In  general,  items  below  the  top  level  need  to  be  changed  only 
rarely  for  a  given  design. 


%  NOTE 
% 

% 

% 

% 

% 

Ts=0. 000005 


17  - 

18 

19  - 

20  “ 
21 

22  “ 

23  - 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 


%  Simulation  step  time,  seconds 
%  Simulation  end  time,  seconds 


Tend=20; 

% 

Nh=4; 

Ng=2; 

% 

Iship0=2762; 

Rh=0.1; 

% 

%  Flywheel  data  for  single  or  multiple  HESM 


% 

% 

% 

% 

% 

% 

% 

% 

% 

% 


%  Number  m  of  HESMs 
%  Number  of  ship  generators 

%  Estimated  constant  current  from  ship's  power  system 
%  HESM  cable  resistance.  Ohms 

used  in  each  HESMm  block 
(read  block's  input  menu)  and  in  yellow  blocks  within  the 
"Measurement”  block  (white  background) . 

Note:  Enter  actual  values  of  FJm  and  v;0m  for  each  HESMm. 

Note:  Match  the  variable  subscript  to  the  subscript  of 
the  HESM  block  actually  used:  e.g.,  FJl  and  wOl 
will  provide  inputs  for  HESMl,  likewise  FJ2  and  w02 
for  HESM2,  etc. 
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36  % 

37  % 

38  % 

39  % 

40  % 

41  % 

42  % 

43  % 

44  % 

45  -  FJl»=1128/Nh; 

46  “  w01=3600; 

47  -  FJ2=1128/Nh; 

48  -  w02=3600; 

49  -  FJ3=1128/Nh; 

50  -  w03=3600; 

51  -  FJ4=1128/Nh; 

52  -  w04=3600; 

53  % 

54  %  PFN  charging 


55 

% 

56 

% 

57 

% 

58 

% 

59 

% 

60 

% 

61 

% 

62 

% 

63 

% 

64 

% 

65  - 

t0=0; 

o 

II 

o 

66  - 

tl=l; 

Pl=23; 

67  - 

t2=4; 

P2=23; 

68  ^ 

t3=4.25; 

P3“0; 

69  - 

t4=5; 

P4*P0; 

Note:  4  HESM  blocks  have  been  provided  with  the  original 
model.  If  Nh>4  new  HESM  blocks  must  be  generated 
copying  them  from  the  other  HESM  blocks  and  making 
the  appropriate  changes  internal  to  the  new  blocks 
(Read  the  block* s  input  menu).  Also  add  additional 
lines  herein  for  the  FJm  and  v;0m  values  needed. 
"Measurement”  block  also  may  be  updated  accordingly 
if  desired. 

%  Flywheel  1  Inertia,  kg*m'"2 
%  Flywheel  1  initial  speed,  RPM 
%  Flywheel  2  Inertia,  kg*m^2 
%  Flywheel  2  initial  speed,  RPM 
%  Flyv;heel  3  Inertia,  kg*m^2 
%  Flywheel  3  initial  speed,  RPM 
%  Flywheel  4  Inertia,  kg*m'^2 
%  Flywheel  4  initial  speed,  RPM 

cycle:  enter  a  sequence  of  5  monotonically  increasing  times 
to  <  tl  <  t2  <  t3  <  t4 
with  corresponding  5  values  of  required  PFN  powers 
PO  PI  P2  P3  P4 

Within  each  time  interval  the  power  curve  will  be 
generated  by  linear  interpolation,  e.g. 

if  tl<t<t2  then  P(t)  =  P1+ (P2-P1) / (t2-tl) * (t-tl) 

If  fewer  than  5  point  pairs  are  sufficient  simply  repeat 
the  last  pair  identical  to  itself  as  needed:  e.g.  if 
only  four  pairs  are  sufficient  make  t4=t3  and  P4=P3. 
Cycle  repeats  up  to  Tend.  Note:  make  Plast=P0. 


70 

71 

72 


% 

% 


END 
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Model  HESM4_PFN1_R1 


IIOT£1‘ 
f&ndJUf1f« 
Me$Mm_PFHljn 
e  funning  »»muU(ion 


HOTC  ti 

PteM«  cfitef  input  djrt  j  for  lA  bl<i<ck« 
with  fionwhMe  bocktprowtd 
by  dovMo-oikking  on  ttve  bioch, 

FoMow  imttnKtiora  given  in  input  menu, 
AhM  look  for  Midtfiorul  inpute  needed 
wttMn  btockt 

(right-eick  erwi  *loelc  Under  MapT 
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Model  HESM2  PFNl  R1 
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